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Preface 


This  rep<M  is  a  primer  on  the  subject  of  structural  dynamics,  oriented 
toward  the  analysis  of  the  intake  structures  that  fall  under  the  authority  of 
the  U.S.  Army  Corps  of  Engineers.  It  is  a  part  of  the  woric  sponsored  by 
the  Computer-Aided  Structural  Engineering  Program  sponsored  by  the 
Directorate.  Headquarters.  U.S.  Army  Corps  of  Engineers  (HQUSACE) 
under  the  Structural  Engineering  Research  Program.  Funds  for  publica¬ 
tion  of  the  report  were  provided  from  those  available  for  the  Computer- 
Aided  Structural  Engineering  (CASE)  project  managed  by  the  Information 
Technology  Laboratory  (ITL)  at  the  U.S.  Army  Corps  of  Engineers  Water¬ 
ways  Experiment  Station  (WES)  at  Vicksburg,  MS.  H.  Wayne  Jones  is  the 
CASE  project  manager.  Technical  monitor  for  the  project  is  Mr.  Lucien 
Guthrie,  HQUSACE. 

The  work  was  performed  at  WES  under  the  direction  of  Dr.  Robert  M. 
Ebeling,  Interdisciplinary  Research  Group,  Computer-Aided  Engineering 
Division  (CAED),  ITL.  The  report  was  prepared  by  Professor  Samuel  E. 
French  of  the  University  of  Tennessee  at  Martin.  Dr.  Ebeling,  and 
Mr.  Ralph  Strom,  North  Pacific  Division,  U.S.  Army  Corps  of  Engineers. 
This  report  is  part  of  a  general  development  of  design  manuals  intended 
for  use  primarily  by  district  personnel  whose  background  does  not  include 
advanc^  dynamics.  The  work  was  accomplished  under  the  general  direc¬ 
tion  of  Dr.  Reed  L.  Mosher,  Acting  Chief,  CAED,  and  the  general  supervi¬ 
sion  of  Dr.  N.  Radhakrishnan,  Director,  ITL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  amUHtt  efihls  report  are  not  to  be  used  for  advenlslng,  pubUeation, 
or  promodoaal  purposes.  Citatloa  ef  trade  names  does  not  coastltuie  au 
offM^  endorsement  or  approval  die  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI  to  SI 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
{Systeme  Internationale)  units  by  applying  the  following  factors: 


HuWiily 

By 

To  Obtain 

angulv  deflraM 

0.01745329 

radians 

cycles  per  second 

1.000 

Hertz 

cycles  per  second 

6.28318531 

radians  per  second 

teet 

0.3048 

meters 

feet  per  second  per  second 

30.4838 

centimeters  per  second  per  second 

gravitational  acceleration 

980.665 

centimeters  per  second  per  second 

gravitational  acceleration 

32.174 

feet  per  second  per  second 

gravitational  aoceleralion 

386.086 

inches  per  second  per  second 

hiches 

2.54 

centimeters 

pounds 

4.4822 

Newtons 

tons 

8.896 

MIonewtons 
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1  Single-Degree-of-Freedom 
Systems 


Introduction 


This  fcport  presents  the  fundamental  concepts  of  the  dynamics  of  in¬ 
take  towers.  The  tfaeoiy  is  quite  general,  however,  and  will  apply  to  multiple- 
degree-of-freedom  (MDOF)  systems  other  than  intake  towers.  The 
ultimate  purpose  is  to  present  the  computer-aided  analysis  of  distributed- 
mass  towers  such  as  the  step-tapered  tower  shown  in  Figure  li.  The  devel¬ 
opment  of  the  theory  therefore  becomes  slanted  toward  distributed-mass 
towers  as  it  progresses.  In  this  chapter  of  the  report  systems  having  a 
single  degree  of  freedom  (SDOF  systems)  are  discussed.  The  single 
lumped-mass  systems  depicted  in  Figure  1  are  typical  examples  of  such 
SDOF  systems. 

The  primary  source  of  dynamic  excitation  of  structures  is  earthquake 
motion.  A  brief  summary  of  earthquake  loading  as  it  is  characterized  in 
terms  of  reqwnse  spectra  is  presented  in  Chapter  2  of  this  report.  Distributed- 
mass  systems  such  as  towers  are  presented  in  Chapter  3  as  an  application 
of  the  Rayleigh  method.  While  the  distributed-mass  systems  of  Chapter  3 
ate  recognized  as  MDOF  systems,  the  Rayleigh  method  allows  them  to  be 
examined  one  degree  of  freedom  at  a  time.  Multistory  lumped-mass  sys¬ 
tems  are  presented  in  Chapter  4.  Such  systems  are  ^OOF  systems  that 
can  also  be  examined  one  degree  of  freedom  at  a  time. 

The  lumped-mass  systems  (but  not  the  distributed-mass  systems) 
shown  in  Figure  1  can  oscillate  in  only  one  way,  that  is,  the  motion  of  the 
single  mass  can  be  only  a  back-and-forth  motion  across  its  initial  at-rest 
position.  Figure  2  is  a  sketch  of  such  an  oscillating  system.  The  position 
of  the  mass  at  any  point  in  time  in  the  system  is  entirely  predictable. 
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Figure  1 .  Common  types  of  oscillating  systems 


Figure  1  shows  several  common  types  of  single-  and  distributed-mass 
oscillating  systems.  If  any  one  of  these  systems  is  given  an  initial  dis¬ 
placement  of  amplitude  uq^  and  then  released,  it  will  oscillate  at  its 
natural  frequency,  continuing  to  oscillate  with  ever-decreasing  amplitude 
(and  ever-decreasing  energy)  until  it  Anally  returns  to  its  at-rest  state. 

The  causes  of  the  decrease  in  amplitude  might  be  such  things  as  hyster¬ 
esis  losses,  inelastic  behavior,  or  any  one  of  many  other  sources  of  energy 
dissipation.  The  cause  could  also  include  man-made  shock  absorbers  de¬ 
liberately  placed  to  reduce  oscillations.  Taken  together,  these  losses  in 


For  convenience,  tymboU  and  abbreviations  are  listed  in  the  notation  (Appendix  B). 
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aaecgy  «e  called  ioMiptef 
widiovt  du^iog, 
die  mau  could  theoreti¬ 
cally  oscillate  forever 
with  iM>  lo^s  in  amplitude. 

Shown  scaematkally, 
damping  effects  are  usu- 
auy  lumped  together  and 
shown  as  a  viscous  shock 
absorber  called  a  dashpot, 
such  as  that  shown  in 
figure  le. 

Figure  3  shows  the  dis¬ 
placement  of  an  oscillat¬ 
ing  mass  both  with  and 
without  damping.  The 
damping  will  affect  only 
the  force  resulting  from 
the  velocity  of  the  mass  in 
motion.  Danqiing  effects 
are  greatest  when  the 
mass  is  at  its  highest  ve¬ 
locity.  As  velocities  de¬ 
crease,  the  damping  also 
decreases,  resulting  in  an 
asymptotic  decrease  in  am- 
plituite  as  shown. 

Figure  3  also  indicates 
the  period  of  oscillation, 

T.  The  period  is  the  time 
in  seconds  for  the  mass  to  complete  one  full  oscillation,  or  cycle.  TYit  fre¬ 
quency  of  oscillation  f  is  another  commonly  used  property;  the  frequency 
is  the  number  of  cycles  completed  in  one  second  (Hz).  Mathematically, 
frequency  is  the  reciprocal  of  the  period  T.  The  rotational  speed  a>  is  an¬ 
other  commonly  used  property  and  is  the  angular  velocity  in  radians  per 
second  where  there  are  2ic  radians  in  one  complete  cycle. 

1  (T 

/  in  cycles/second  =*  := - 7  ' 

^  '  T  seconds 


Figure  2.  Oscillations  of  a  lumped-mass 
system 


(D  in  radians/ second  = 


2ic  radians/c 


T  seconds 


(2) 


Figure  3  shows  that  the  period  of  oscillation  T  for  this  structure  is  es¬ 
sentially  the  same  for  both  damped  and  undamped  oscillations.  At  very 
high  levels  of  damping,  there  will  in  fact  be  a  slight  elongation  in  the 
period  T  (Toy^PO)  >  T^i^dampeo)'  f®*”  ordinary  structures  working 
at  elastic  levels  of  stress,  the  level  of  damping  is  so  low  that  any  such 
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Figure  3.  Damped  oscillations 


difference  can  be  safely  ignored.  The  dominant  effect  of  damping  is  to  re¬ 
duce  the  velocity  of  oscillation,  t'^*reby  producing  the  asymptotic  reduc¬ 
tion  in  amplitude  shown  in  Figure  3. 

The  angular  velocity  m  does  not  always  appear  physically  in  an  oscillat¬ 
ing  system,  but  the  symbol  m  usually  does.  Figure  4,  in  which  an  imagin¬ 
ary  wheel  is  added  to  a  laterally  oscillating  system,  shows  a  means  to 
visual!.  '  a  rotational  velocity  when  none  actually  exists.  The  rotational 
velocity  of  the  imaginary  wheel  is  o). 


Figure  4.  Rotational  velocity  for  a  linear  oscillation 
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Equation  of  Motion 


The  equation  of  motion  for  lumped-mass  systems  having  a  single  de¬ 
gree  of  freedom  (termed  SDOF  systems)  is  derived  from  Newton’s  second 
law,  force  =  mass  times  acceleration.  Consider  the  system  shown  in  Fig¬ 
ure  5  in  which  the  oscillating  mass  has  been  removed  as  a  free  body  at 
some  instant  in  time  t.  An  externally  applied  dynamic  force  p(t)  is  acting 
on  the  system.  The  force  p(t)  varies  with  time  and  is  called  a.  forcing 
function. 


Figure  5.  Common  SDOF  oscillating  system 


The  free  body  shown  in  Figure  5  is  subject  to  four  forces,  each  one 
varying  with  time. 

a.  The  inertial  force  (f^  =  mil)  is  given  by  D’Alembert’s  principle  to  be 
a  force  acting  through  the  center  of  mass,  opposite  in  direction  to 
the  imposed  motion. 

b.  The  restoring  force  (f|^  =  ku)  is  the  stiffness  k  (in  pounds  of  force 
per  unit  of  displacement)  times  the  displacement  u. 

c.  The  damping  force  (f^  =  cu)  is  the  damping  constant  c  (in  pounds  of 
foree  per  unit  of  velocity)  times  the  velocity  u  at  any  time  t. 

d.  The  fourth  force  is  that  of  the  forcing  function  [p(t)],  which  may  or 
may  not  exist  at  any  given  time. 

The  sum  of  forces  yields  the  equation  of  basic  force  equilibrium: 


-/i  -fk  -fd  +  P(#)  =  0 


(3) 


CiMptarl  8tngl>-Degr—  of-Fr— dom  Systenw 


AU  of  dme  foices  may  be  expiested  in  terms  of  the  displacement  u, 
whme,  as  noted  earlier. 


//•«*» 

(4a) 

/^»  ci 

(4b) 

4 -to 

(4c) 

Substicnting  Equations  4a.  b.  and  c  into  Equation  3  yields  the  basic 
equation  of  motion  for  an  SDOF  system  subjected  to  an  externally  aiqtlied 
ftncing  function: 

mil*  +  cii  -f  s  p(t)  (S) 

This  same  equation  applies  equally  to  all  of  the  SDOF  systems  shown 
in  Hgure  1.  Fm  the  rotational  oscill^tms  of  Figure  Ih,  rotational  diS' 
placement  n  would  have  to  be  substituted  for  lateral  displacement  u. 

In  all  itf  the  liuiq)ed-mass  systems  of  figure  1.  the  mass  of  the  system 
is  understood  to  be  rigid.  Ccwskter,  for  example,  the  idealized  one-story, 
single-bay  building  fiame  of  Hgtite  Id.  Under  the  indicated  displace¬ 
ment.  it  is  assumed  that  the  rigid  mass  does  not  undergo  any  bending:  only 
the  columns  will  experience  bending.  Therefixe,  all  of  the  restoring  force 
f|.  in  this  system  must  come  from  bending  in  the  coluirms. 

The  equation  of  motion  (Equaticm  5)  for  the  systems  of  Figure  1  in¬ 
cludes  a  forcing  fimctiim  p(t)  applied  to  the  mass  m.  An  earthquake,  how¬ 
ever,  does  not  produce  such  a  tangible  force  acting  on  the  system.  Rather, 
an  earthquake  will  produce  displacements  of  the  foundations  of  the  struc¬ 
ture.  It  is  Aese  diqilacenients,  rather  than  uagible  forces,  which  the 
earth^iake  inqwses  on  the  system. 

The  results  of  such  earthquake-induced  displacements  are  shown  in  the 
free  body  of  Figure  6.  The  earthquake  ground  motirms  ate  shown  as  Ug(t), 
which  ate  random  motioiu  in  time.  The  total  displacement  of  the  mass  rel¬ 
ative  to  its  at-test  position,  U|(t),  is  shown  in  Figure  6a.  The  displacement 
u  again  represents  the  displacement  of  the  structural  mass  relative  to  the 
base  of  the  structure  at  any  time  t.  The  geometry  of  Hgure  6a  shows  that 
the  total  diqrlacement  U|(t)  is  given  by 

«i(t)  *  itg(0  +  "(f)  (6) 

Figure  6b  shows  the  forces  associated  with  the  free  body  of  Figure  6a. 
fa  this  case,  however,  the  inertial  force  f)  created  by  tlm  total  acceleration 
Ot(t)  becomes: 

/f  m  nmfi)  *  mii^r)  +  inii  (7) 


e 
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Figure  6.  Relative  displacements  of  SDOF  systems 


Forces  are  summed  as  before,  yielding  the  equation  of  force  equilibrium: 

where 

fd  ®cu 
f]^  =  ku 

These  values  of  displacements  are  substimted  into  Equation  8,  yielding: 
ma^(r)  +  mil  +  cii  +  ku  =  0  (9) 

Rearranging  terms  produces  the  final  form  of  the  equation  of  motion  for  a 
structure  subjected  to  earthquake  displacements  at  its  base: 

nui  +  cu  +  ku  -m«^(r)  (10) 

Equation  10  is  essentially  the  same  as  Equation  5.  with  the  forcing  func- 
tiiHi  being  the  effective  force  -m(l  (t).  The  amount  of  force  exerted  on  a 
structure  is  therefore  a  function  of  the  mass  of  the  structure  itself.  Light¬ 
weight  steel  or  timber  structures  will  experience  much  less  force  from 
earthquake  motions  than  will  heavier  concrete  or  masonry  structures. 

Equation  10  is  the  general  equation  of  motion  for  a  single-mass  struc¬ 
ture  subjected  to  a  dynamic  forcing  function  such  as  an  earthquake.  A  typ¬ 
ical  exanqile  of  such  a  structure  is  shown  schematically  in  Figure  7. 

Bgure  7  shows  that  under  the  actual  excitation,  the  earth  is  moving  relative 
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Figure  7.  SDOF  system  under  earthquake  forcing  function 


to  the  mass  whereas  under  the  equivalent  excitation,  the  mass  is  moving 
relative  to  the  earth. 

The  major  part  of  this  report  is  devoted  to  the  solution  of  Equation  10. 
As  with  any  differential  equation,  the  solution  of  Equation  10  consists  of 
two  parts: 

a.  The  general  solution  of  the  homogeneous  equation,  that  is,  a 
solution  with  the  forcing  function  0^(0  s  0  and  the  structure 
undergoing  unforced  (free)  oscillations. 

b.  A  particular  solution  of  the  general  equation,  that  is,  a  solution  with 
the  forcing  function  in  place  and  the  structure  being  forced  to  undergo 
the  same  displacements  as  the  forcing  fiinctitm,  but  with  a  time  lag. 

The  solution  for  the  unforced  free  oscillations  will  be  developed  first. 


Undamped  Free  Oscillations 

The  equation  of  motion  given  by  the  homogeneous  portion  of  Equa- 
timi  10  will  be  solved  flrst  without  any  damping.  Damping  will  then  be 
adcted,  and  the  final  form  of  the  solution  will  be  developed. 

At  some  instant  t  s  tQ,  an  SDOF  system  is  assumed  to  have  a  known  dis¬ 
placement  Uq  and  a  known  velocity  .  Thereafter,  the  system  is  allowed  to 
oscillate  freely  without  any  further  fwcing  function  being  applied.  The  equa¬ 
tion  of  motion  for  the  undunped  free  oscillations  is  then  given  by: 

mil  4-  ku  =  0  (11a) 


or, 

ii  +  (k/m)u  s  0 


(lib) 
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Tie  totatkNi  of  oi^aasy  ^ffuential  eqn^tm  is: 


M  «  A  8m(Qor)  +  B  cosCfor)  (12) 

where  <d  is  dM  c<mttant  circular  frequency  of  the  solution  in  radians  per 
second  and  A  and  B  are  the  constants  of  integration. 

The  circular  frequency  o>  is  found  from  the  solution  as: 

<0  »  VtTiir  (13) 

The  natural  period  of  oscillation  T  can  now  be  determined,  even  before 
the  constants  of  integration  have  been  evaluated: 

r  =  2it  -^mTk  (14) 

Equation  14  shows  that  the  natural  period  of  oscillation  of  an  un- 
danqMsd  SDOF  system  is  dependent  only  on  the  mass  of  the  system  and 
the  stiffness  of  the  system.  No  boundary  conditions  or  any  other  factors 
affect  the  undamped  natural  period  T. 

The  constants  of  integration  A  and  B  can  be  found  as  usual  by  evaluat¬ 
ing  the  boundary  conditions: 

At  time  t  =  0,  u  =  Uq  =  A  sin[(D(0)]  +  B  cos[o>(0)] 

«0  =  (15) 

The  velocity  it  at  any  time  t  is  given  by: 

u  =  A  00  cos  (oor)  -  B  00  sin  (oor) 

At  time  r  =  0.  «  =  =  Aoo 

Uq/oo  =  A  (16) 

The  final  solution  is  then  given  by: 

«  s  11^  cos  (oor)  +  (miq/oo)  sin  (oor)  (17) 

The  maximum  amplitude  given  by  Equation  17  will  occur  when  the  deriva¬ 
tive  of  die  ampliti^  (the  velocity)  is  zero: 

H  s  0  =  sin  (oor)  +  (hq/oo)  oocos  (oor) 
or. 


I  0  s  -  ii^  sin  (oor)  -t-  («g/oo)  cos  (oor)  (18) 

I 

i 
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DampMi  Frao  Oscillations 


The  foregoing  solution  for  ftee  oscillations  did  not  include  damping. 
When  damping  is  included  (still  without  a  forcing  function),  the  homoge¬ 
neous  part  of  Equation  10  becomes: 

mil  -f  cu  +  Au  s  0  (22) 

The  equation  is  divided  by  m  to  put  it  in  standard  form: 

u  +  2p()iM  +  a?u  =  0  (23a) 

where  the  constant 

P  =  c/2nno  (23b) 


For  p  <  1,  the  solution  of  Equation  23a  is: 


ii  =  e-P“ 


Uq  cos  (UqT  + 


it  + 


sin 


(24) 


The  motion  described  by  Equation  24  for  P  <  1  is  summarized  in  Figure  9. 
It  is  a  periodic  function  having  an  exponentially  decaying  amplitude  and  a 
damped  circular  frequency  (Oq: 


<D|,  =  m  Vl  -  p2  =  V;fc(l  --  pV« 


(25) 


where  co  is.  as  usual,  the  circular  frequency  of  the  undamped  system. 

The  damped  period  for  the  motion  of  Equation  24  is  found  from  Equa¬ 
tion  25: 

Tjj  =  tM\  -  p2  (26) 


where  T  is.  as  previously  deEned.  the  natural  period  of  oscillation  of  the 
undamped  system. 

If  the  factor  P  in  Equation  26  could  be  made  equal  to  1.  the  natural  pe¬ 
riod  of  the  danced  system  would  become  infinite.  Such  a  system  would 
simply  return  to  its  original  at-rest  position  following  a  displacement,  i.e., 
there  would  be  no  oscillations.  For  the  case  in  which  P  =  1 ,  the  damping 
coefficient  c  as  given  by  Equation  23  is  its  value  of  critical  damping: 

=  ^crttical=  2m®  (27) 


Chapter  1  Singte-Dagraa-«f-Fraadom  Systems 


Further  examination  of  Equation  23a  reveals  that  the  factor  P  may  be 
regarded  as  the  fraction  of  critical  damping  in  a  system,  i.e.: 

P  =  c/2ma)  =  c/c^rttic<a  (28) 

For  typical  structures  woilting  in  their  elastic  ranges,  the  factor  P  has 
been  found  to  be  less  than  0.1,  or,  stated  another  way,  the  maximum 
amount  of  damping  in  a  typical  structure  will  normally  be  less  than  10  per¬ 
cent  of  its  critical  danq>ing.  Substitution  of  this  maximum  value  of  P  =  0. 1 
into  Equation  26  indicates  that  for  typical  structures,  the  effects  of  damping 
on  the  period  of  oscillation  are  indeed  negligible  and  can  be  ignored,  i.e.: 
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r-r^ 


(29) 


Forced  Oscillations  with  a  Dynamic  Forcing 
Function 


As  stated  earlier,  the  intended  solution  of  Equation  10  would  consist  of 
the  general  solution  of  the  homogeneous  equation,  both  with  and  without 
damping,  plus  a  particular  solution  of  the  general  equation,  both  with  and 
without  damping.  Equations  17  and  24  satisfy  the  first  part  of  this  solu¬ 
tion.  Hiey  are  the  required  general  solutions  of  the  homogeneous  equa- 
ti<m.  It  remains  now  to  And  a  particular  solution  of  Equation  10,  both 
with  and  without  damping,  when  the  system  is  subjected  to  a  dynamic 
forcing  function. 

When  a  forcing  function  is  included.  Equation  10  takes  its  most  general 
form: 

+  CM  -t-  ku  s  -niujjt)  (10) 

The  forcing  function  mfl  (t)  is  the  mass  of  the  structure  times  the  ground 
acceleration:  physically,  the  result  is  a  force.  The  nature  of  this  force  has 
yet  to  be  deA^. 

Two  cases  of  dynamic  forcing  functions  will  be  considered  for  examina¬ 
tion.  one  being  a  harmonic  variation  and  the  other  being  the  random  varia¬ 
tion  produced  by  an  earthquake.  The  first  case  considers  a  forcing 
function  that  has  a  harmonic  variation  similar  to  that  of  the  natural  oscilla¬ 
tions  of  the  structure.  The  forcing  function,  however,  is  assumed  to  have 
some  frequency  which  is  not  necessarily  the  same  as  the  natural  fre¬ 
quency  of  the  structure: 

-mu  as  Pq  sin  mr  (30) 

where  Pq  is  the  maximum  magnitude  of  the  forcing  function. 

When  a  harmonic  forcing  function  such  as  that  of  Equation  30  is  ap¬ 
plied  to  a  damped  SDOF  structure,  the  response  of  the  structure  will  con¬ 
sist  of  two  parts,  as  noted  earlier.  The  Arst  part,  the  free  oscillation  part, 
is  that  given  by  Equation  24.  This  free  oscillation  soon  decays,  however, 
leaving  only  the  second  part.  The  second  part  is  the  continuous  steady- 
state  response  that  is  being  imposed  on  the  system  by  the  forcing  function. 
This  ste^-state  response  is  the  particular  solution  of  Equation  10.  The 
forcing  function  is  physically  forcing  the  SDOF  system  to  follow  the  mo¬ 
tions  of  the  forcing  function. 
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However,  the  forcing  function  is  not  acting  on  the  oscillating  mass. 
The  forcing  function  acts  at  the  base  of  the  columns  as  shown  earlier  in 
Figure  7.  The  motion  of  the  mass  is  being  created  by  a  back-and-forth 
ground  motion  at  the  base  of  the  columns.  At  any  time  t,  the  force  im¬ 
parted  to  the  columns  is  that  given  by  Equation  30. 

The  displacements  of  the  mass  will  follow  the  same  harmonic  displace 
ments  as  the  ground  motion.  There  will  be  a  lag  in  time,  however,  be¬ 
tween  the  motion  of  the  ground  and  the  motion  of  the  nuiss.  The  lag 
occurs  as  the  columns  develop  enough  lateral  force  to  produce  (or  to 
change)  the  motion  in  the  mass.  It  is  well  to  note  also  that  the  amplitude 
of  motion  of  the  mass  need  not  be  equal  to  the  amplitude  of  motion  of  the 
ground.  Such  a  difference  is  indicated  in  Figure  10. 


The  displacement  u  of  the  oscillating  mass  at  any  time  t  is; 

“  =  “msx  (31) 

where  9  is  the  lagging  phase  angle  between  the  motion  of  the  ground  and 
the  motion  of  the  mass. 

The  solution  for  the  response  of  the  system  to  various  excitation  fre¬ 
quencies  m  is  mathematically  quite  complex,  far  beyond  the  scope  of  an 
elementary  treatment  such  as  this  (see  Paz,  Structural  Dynamics  1991). 

The  end  result  of  the  solution,  however,  is  relatively  simple.  The  displace¬ 
ment  u  at  any  time  t  is  given  by  the  solution  as: 

“  =  “nuu  ~  (32a) 
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(32b) 


•‘ti  Vn  -  +  [2fi(&/a)f 

aod  the  lagging  phase  angle  6  is  given  by 

e  .  (32c) 

1  -  (&/to)^ 

In  EqnatitMis  32a  and  32b,  the  displacement  is  the  displacement  that 
the  maximum  force  p^  would  produce  if  the  force  were  to  be  applied  as  a 
static  force  (i.e.,  u,^ »  P</h).  All  other  symbols  in  Equations  32a,  b,  and  c 
were  introduced  and  denned  earlier. 

Equation  32b  reveals  that  there  are  two  primary  variables.  The  first 
variable  is  Um«y/u,^,  which  indicates  the  relative  level  of  magnification 
of  displacement  at  various  forcing  frequencies  &.  The  second  variable  is 
S/to ,  which  is  the  ratio  of  the  forcing  frequency  S  to  the  natural  fre¬ 
quency  (o.  (The  third  variable  fi,  the  ratio  of  critical  damping,  can  be  as¬ 
signed  whatever  value  the  physical  conditions  warrant.) 

Figure  1 1  inesents  a  graph  of  Equation  32b  for  various  levels  of  damp¬ 
ing.  The  most  significant  feature  of  the  graph  is,  of  course,  the  sharp  in¬ 
crease  in  displacement  as  the  frequency  of  the  forcing  function 
apinoaches  the  natural  frequency  of  the  system.  A  second  significant  fea¬ 
ture  of  the  gt^h  is  the  effect  of  the  damping.  A  large  amount  of  damping 
is  required  if  a  pronounced  reduction  in  the  peaking  effect  is  required. 

Hgure  1 1  reveals  the  phenomenon  of  resonance.  The  peak  value  of  dis¬ 
placement  given  by  Equation  32b  occurs  when  the  forcing  frequency  S  is: 

In  comparison,  the  damped  natural  frequency  a  is: 

csp  =  mVl  -  ^  (34) 

The  resonant  frequency  is  therefore  not  exactly  equal  to  the  damped 
natural  frequency  of  the  system,  but  for  low  values  of  damping,  the 
difference  is  quite  small.  In  general,  this  small  difference  is  ignored, 
and  the  resonant  frequency  is  taken  to  be  essentially  equal  to  the  natu¬ 
ral  frequency. 

Eisofar  as  the  structural  design  is  concerned,  it  is  considered  prudent  to 
keep  the  natural  frequency  of  a  structure  as  far  away  as  possible  from  the 
frequency  m  of  any  known  dynamic  loading.  Since  the  natural  frequency 
of  the  structure  is  a  frinction  only  of  stiffness  and  mass: 
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For  the  general  case  of  loading  in  which  the  forcing  function  is  that  of 
random  earthquake  motions,  the  equation  of  motion.  Equation  10,  remains 
in  its  most  geimral  form: 


mm  +  at  +  ku  -  -mm^(t)  (10) 

In  its  standard  integrable  form.  Equation  10  becomes: 
ii  +  2PaM  +  cD?ii  s* 

The  problem  of  a  damped  SDOF  system  shaken  by  a  time-varying  ground 
acceleration  is  equivalent  to  the  problem  of  a  damped  SDOF  system  rest¬ 
ing  on  a  fixed  baM  and  being  subjected  to  a  time-dependent  force  p(t)  of 
magnitude  -m  •  ii^(r),  as  shown  in  Figure  12. 


The  total  dynamic  response  of  the  SDOF  system  problem  is  computed 
in  two  steps: 

Step  1)  Solve  for  the  relative  response  of  the  damped  SDOF  system 
as  governed  by  the  ordinary  differential  equation. 

Equation  35. 
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Step  2)  Sum  the  relative  response  with  the  motion  of  the  ground  to 
obtain  the  total  response. 

Step  1  assumes  that  the  contribution  of  the  general  solution  is  negligible 
compared  with  the  contribution  of  the  particular  solution. 

Few  simple  harmonic  ground  accelerations,  e.g.,  (ug  =  constant  ■  sin  (a>/)], 
closed-form  solutions  to  Equation  35  are  available  in  numerous  textbooks 
on  both  mechanical  vibrations  and  structural  dynamics.  This  procedure  is 
impractical  for  solving  earthquake  engineering  problems  due  to  the  irregu¬ 
lar  nature  of  ground  acceleration/time  histories. 

A  second  procedure  used  to  solve  for  the  relative  displacement  of  the 
SDOF  system  involves  the  representation  of  the  load/time  history 
p(t)  =  -m  ■  iigU)  as  a  series  of  impulse  loadings  P(t)  applied  to  the  SDOF 
system  for  infinitesimal  time  intervals  dx.  Duhamel’s  integral  for  a 
damped  SDOF  system  is  given  by: 


«(/)  =  -(1/(0^  f  Ug(x)e-^*-'^>  sin  [<o^(t  -  x)]dt 


(36) 


where 


ODD  s  the  damped  angular  frequency  of  vibration. 

<D  s  the  undamped  angular  frequency  of  vibration. 

P  =  the  fraction  of  critical  damping. 

The  irregular  forms  of  acceleration/tinK  histories  require  numerical  solu¬ 
tions  to  be  used  to  solve  Duhamel’s  integral  (e.g.,  Paz  1991,  or  Clough 
and  Penzien  1$193). 

In  usual  applications  to  earthquake  engineering  problems,  numerical 
methods  are  used  to  solve  Equations  35  and  36  for  the  relative  displace¬ 
ment  of  the  SDOF  mass  because  of  the  irregular  nature  of  ground  accelera¬ 
tion/time  histories.  In  general,  there  are  two  categories  of  numerical 
methods  used  for  solving  the  dynamic  equilibrium  equation:  1 )  direct  inte¬ 
gration  methods,  and  2)  frequency  domain  methods.  The  reader  is  re¬ 
ferred  to  books  on  structural  dynamics,  e.g.,  Paz  1991,  Clough  and 
Penzien  1993,  or  Ebeling  1992  for  a  description  of  these  two  methods. 

The  solution  of  Duhamel’s  equation  with  p  =  0  yields  a  useful  relation¬ 
ship  between  displacement,  velocity,  and  acceleration  (Appendix  A,  Ebel¬ 
ing  1992). 

it  =  (ou  (37a) 
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M  K  coi^ 


(37b) 


It  is  important  to  recognize  that  the  relationships  given  by  Equations  37a 
and  b  are  exact  only  when  damping  is  zero.  The  velocity  u  computed 
from  Equation  37a  is  called  pseudovelocity  and  the  acceleration  Q  com¬ 
puted  Arom  Equation  37b  is  called  pseudoacceleration.  Equations  37a 
and  b  define  the  interrelationship  tetween  spectral  displacement  u, 
pseudovelocity  u,  and  pseudoacceleration  Q,  three  terms  to  be  discussed  in 
detail  in  Chapter  2. 
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2  Design  Response  Spectra 


Introduction 


Earthquake  loading  is  typically  represented  in  a  dynamic  analysis 
either  by  a  ground  acceleration/time  history  or  by  a  response  spectrum.  A 
response  spectrum  is  a  graphical  relationship  of  maximum  values  of  accel¬ 
eration.  velocity,  and  displacement  response  of  SDOF  systems  having  a 
natural  frequency  f.  The  graph  is  drawn  over  the  usual  range  of  frequen¬ 
cies  for  elastic  SDOF  systems. 

To  prepare  such  a  response  spectrum.  Equation  35  is  evaluated  for  the 
relative  accelerations  ii  of  an  SDOF  system  during  a  particular  earthquake 
excitation.  The  evaluation  is  made  for  the  frequency  f  (or  period  T)  of  the 
SDOF  system  and  its  level  of  damping  p.  As  indicated  in  Figure  13.  the 
absolute  acceleration  response  spectrum  SA  is  then  computed  as  the  maxi¬ 
mum  absolute  value  of  the  sum  of  this  computed  relative  acceleration/time 
history  for  the  SDOF  system  plus  the  ground  acceleration/time  history 
from  step  2  of  Figure  1 2. 

Also  indicated  in  Figure  13,  the  solution  is  then  searched  over  the  en¬ 
tire  time  history  for  the  maximum  absolute  acceleration  SA  of  the  SDOF 
system,  which  is  found  and  recorded.  The  corresponding  value  of 
pseudoacceleration  S^  is  also  computed  and  recorded,  where  S^  = 

The  procedure  is  performed  repeatedly  at  close  intervals  of  frequency. 

The  resulting  sets  of  values  for  frequency  f  and  pseudoacceleration  S^^  are 
then  plotted  as  the  spectra  of  responses  for  the  various  SDOF  systems.  A 
key  feature  of  a  response  spectrum  is  that  it  has  been  made  independent  of 
time. 

As  more  recordings  of  earthquake  acceleration/time  histories  became 
available  during  the  1970’s,  statistical  analyses  of  response  spectra  of 
earthquake  ground  motions  were  conducted  by  several  groups  of  earth¬ 
quake  engineers  and  seismologists.  The  results  of  some  of  the  early 
ground  motion  studies  are  summarized  in  Seed,  Ugas,  and  Lysmer  (1976), 
Mohraz  (1976),  and  Newmark  and  Hall  (1982).  Their  studies  showed  that 
the  spectral  frequency  contents  of  the  recorded  accelerograms  were  depen¬ 
dent  on  the  earthquake  magnitude,  distance  from  causative  fault  to  site. 
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Figure  14.  Smoothing  of  raw  response  data 


the  peak  values.  The  overall  accuracy  increases,  of  course,  as  the  number 
of  records  increases. 

With  an  adequate  number  of  records  (see  the  appendix  of  EM  1 1 10-2- 
6050),  such  smooth  curves  can  be  developed  well  within  the  range  of  accu 
racy  of  the  earthquake  projections  themselves.  All  of  the  smooth-shaped 
spectra  are  drawn  using  this  technique. 


Smooth-Shaped,  Broad-Band  Design  Besponse 
Spectra 

Early  ground  motion  studies  identified  the  factors  affecting  the  re¬ 
sponse  spectra  for  earthquake  motions,  characterizing  shape  or  the  fre¬ 
quency  content  of  the  earthquake  spectrum  for  the  category  under  study 
and  developing  smooth,  broad-band  spectra  for  use  in  the  design  of  struc¬ 
tures  for  earthquake  loadings.  A  broad-band  spectrum  ensures  that  suffi¬ 
cient  seismic  energy  is  delivered  to  all  frequencies. 
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•  Th«  Newmufe  nd  Hall  vespoate  ipeetra. 

•  The  Apidied  Technology  Council  ATC  3-06  response  spectra. 

BoA  spectra  era  nonsite-specific  and  both  can  be  an>lied  in  general  strnc 
tnral  iqiplications. 


Newmark  and  Hall  Responae  Spectra 

In  the  Newnurk  and  Hall  approach  for  developing  design  response 
spectra,  it  is  necessary  first  to  identify  the  peak  ground  acceleration 
(POA),  peak  ground  velocity  (POV).  and  ground  displacement 
(POD)  for  the  design  earthquidce.  These  peak  values  are  then  multiplied 
by  Impropriate  spectral  ampliflcation  factors  to  obtain  the  corresponding 
diqrlacements,  velocities,  and  accelerations  of  SDOP  systems  throughout 
a  range  of  natural  frequencies.  The  results  are  presented  graphically  in  a 
particular  type  of  griqih  called  a  tripartite  ^aph.  Figure  15  is  an  example 
of  such  a  design  response  spectrum  adiq)^  f'ro"*  Newmaik  and  Hall  1982. 


Figure  15.  Newmark  and  Hall  design  response  spectrum  for  a  site  having 
a  stiff  soil 


There  are  four  logarithmic  scales  used  in  Figure  IS  to  show  the  follow¬ 
ing  three  quantities  as  functions  of  the  frequency  f. 

•  D'  is  the  maximum  relative  displacenwnt  between  the  ground  and 
the  mass  of  an  SDOF  system  having  a  natural  frequency  f,  at 
whatever  point  in  time  it  occurs. 

•  V'  is  the  maximum  pseudovelocity  cs^D'  of  the  SDOF  mass. 

•  A'  is  the  maximum  pseudoacceleration  of  the  SDOF  mass. 

Since  all  four  quantities  are  interrelated,  the  tripartite  plot  provides  a  con¬ 
venient  means  to  show  the  entire  interrelationship  in  a  single  curve.  It  can 
be  shown  (Hudson  1979  or  Ebeling  1992)  that  for  low  levels  of  damping, 
the  pseudovelocity  and  the  pseudoacceleration  are  nearly  equal  to  the  max¬ 
imum  relative  velocity  and  maximum  relative  acceleration.  This  approxi¬ 
mation  applies  over  most  of  the  usual  frequency  range,  with  the 
pseudoacceleration  being  more  accurate  than  the  pseudovelocity.  (For  ad¬ 
ditional  details  regarding  these  issues,  see  Ebeling  1992.) 

For  reference,  the  maximum  ground  motion  is  shown  in  dashed  lines 
on  the  graph  of  Figure  1 S.  The  ground  motion  shows  only  values  of  PGD, 
PGV,  and  PGA  without  regard  to  time. 

In  order  to  draw  a  SDOF  response  curve,  a  ground  motion  curve  must 
be  established,  which  means  that  the  values  of  PGD,  PGV,  and  PGA  must 
be  established.  To  establish  these  peak  values  of  ground  motion,  it  is  nec¬ 
essary  first  to  adopt  a  design  basis  earthquake.  There  are  two  minimum 
design  basis  earthquakes  in  use  by  the  Army  Corps  of  Engineers  for  intake 
towers  (ER  1 1 10-2-1806  and  EM  1 1 10-2-2401): 

•  Operational  Basis  Earthquake  (OBE):  The  OBE  is  the  level  of 
ground  motion  for  which  the  structure  is  expected  to  remain 
functional  with  little  or  no  damage.  Ordinarily,  the  OBE  is  defined 
as  a  ground  motion  having  a  SO  percent  probability  of  exceedance 
during  the  design  life  of  100  years  (a  144-year  return  period).  The 
associated  performance  level  is  normally  the  requirement  that  the 
structure  will  function  within  the  elastic  range  with  little  or  no 
damage  and  without  interruption  of  function.  Because  the  purpose 
of  the  OBE  is  to  protect  against  economic  losses  from  damage  or 
loss  of  service,  alternative  choices  of  return  period  for  the  OBE 
may  be  made  on  the  basis  of  an  economic  analysis.  The  OBE  is 
normally  based  on  a  probabilistic  site  hazard  analysis  (PSHA). 

•  Maximum  Design  Earthquake  (MDE):  The  MDE  is  the  maximum 
level  of  ground  motion  for  which  the  structure  is  designed  or 
evaluated.  The  associated  performance  level  is  the  requirement 
that  the  structure  perform  without  catastrophic  failure,  such  as 
uncontrolled  release  of  a  reservoir,  although  severe  damage  or 
economic  loss  may  be  tolerated.  For  critical  structures  (refer  to 
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ftR  1110>2-1806).  Ike  Mm  U  the  seme  es  the  nuudmun  ciedible 
Midiqiieke  (hTCB).  For  noaeritical  stmctores  (refer  to  ER  1 1 10-2- 
1806).  the  MDE  mey  be  chosen  es  e  lesser  earthqneke  dian  die 
MCE  to  provide  ^nomicel  designs  meeting  appropriate  safeQr 
standards.  For  noncritkal  stnictnres,  the  MDE  can  be  determined 
by  aPSHA. 

The  maximum  credible  earthquake  (MCE)  is  defined  as  the  greatest 
earthquake  that  can  reasonably  be  expected  to  be  generated  by  a  specific 
source,  on  the  basis  of  seismological  and  geological  evidence.  Since  a 
project  site  may  be  affected  by  earthquakes  generated  by  various  sources, 
each  with  its  own  fault  mechanism,  maximum  earthquake  magnitude,  and 
distance  from  the  site,  multiple  MCE*s  may  be  defined  for  the  site,  each 
with  characteristic  ground  motion  parameters  and  spectral  shape.  The 
MCE  is  determined  by  a  Deterministic  Seismic  Hai^  Analysis  (DSHA). 

Moderate  levels  of  damage  may  be  acceptable  in  existing  intake  towers 
if  it  can  be  demonstrated  that  the  structure  will  continue  to  function  fol¬ 
lowing  an  OBE  event  The  acceptance  of  a  higher  level  of  damage  in  ex¬ 
isting  intake  towers  will  depend  on  the  cost  of  retrofitting  to  preclude 
damage  versus  the  cost  of  repairs  following  an  OBE  event. 

There  are  two  special  cases  for  defining  the  MDE  to  be  the  same  as  the 
MCE.  In  the  special  case  in  which  failure  of  the  tower  due  to  an  earth¬ 
quake  can  lead  to  failure  of  the  dam  and  loss  of  the  reservoir,  the  MDE 
should  be  the  maximum  ciedible  earthquake  (MCE). 

The  other  special  case  in  which  the  MDE  is  defined  as  the  MCE  occurs 
when  the  intake  tower  is  required  to  operate  after  a  severe  earthquake. 

The  tower  may  be  damaged  but  its  ability  to  function  is  not  impaired  by 
damage  sustained  during  ground  motions  in  excess  of  the  minimum  design 
basis  earthquake.  The  MDE  is  the  MCE  for  this  case  so  that  the 
postearthquake  functionality  of  the  intake  structure  can  be  demonstrated. 

When  site-specific  information  is  not  available,  Newmark  and  Hall  rec¬ 
ommend  that  a  v/a  ratio  of  48  in/sec/g  (122  cm/sec/g)  be  used  for  a  compe¬ 
tent  (stiff)  soil.  For  rock,  they  recommend  a  v/a  ratio  of  36  in/sec/g 
(91  cm/sw/g).  To  ensure  that  the  spectrum  includes  an  adequate  band 
width  of  frequencies,  they  recommend  that  ad/v^  be  taken  at  6.  (In  these 
recommendations,  a,  v,  and  d  are  the  values  of  the  peak  ground  accelera¬ 
tion,  peak  ground  velocity,  and  peak  ground  displacement,  respectively.) 
From  these  recommended  values,  the  value  of  d/a  is  found  to  be  90  cm/g 
for  competent  soil  and  SI  cm/g  for  rock. 

The  smooth  elastic  design  spectrum  shown  in  Figure  12  is  an  84^  per¬ 
centile  spectrum  developed  for  a  site  located  on  a  competent  soil  with  an 
estimated  PGA  value  of  a  =  O.Sg.  In  the  spectrum  of  Figure  15,  a  value  of 
122  cm/sec/g  is  used  for  v/a,  and  a  value  of  90  is  used  for  d/a.  For  a  value 
of  a  s  O.Sg,  the  corresponding  value  for  velocity  v  is  61  cm/sec  and  for 
displacement  d  is  45  cm. 
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In  an  earlier  paper.  Newmark  and  Hall  (1978)  provided  values  of  ampli- 
Hcation  factors  for  the  different  parts  of  the  spectrum.  These  spectrum 
amplification  factors  are  shown  in  Table  1  for  various  damping  ratios. 

For  the  MDE,  the  spectrum  amplification  factors  at  5  percent  damping  are 
2.71, 2.30,  and  2.01  for  acceleration  A,  velocity  V,  and  displacement  D.  re¬ 
spectively,  of  die  SDOF  system. 


Tabtol 

Spoctnmi  Ampllficafion  Factort  for  Horizontal  Elaatic  ftoaponaa 
(Sourca:  Adi^lad  from  Nawniark  and  Hall  1962) 

OTWWi 

OmalaiM^il^ 

IMIan(80%)  1 

A 

V 

0 

A 

V 

D 

0.8 

8.10 

3.84 

3.04 

3.68 

239 

2.01 

1.0 

4.38 

3.38 

2.73 

3A1 

2A1 

1.82 

2.0 

3.08 

2.02 

2.42 

2.74 

2.03 

1.83 

3.0 

3.24 

2.84 

2.24 

2.48 

1.88 

1.82 

8.0 

2.71 

2.30 

2.01 

2.12 

1.85 

1.39 

7.0 

2.38 

2.06 

IAS 

1.89 

1.81 

1.29 

10.0 

1.98 

1.84 

1.69 

1.64 

137 

1.20 

20.0 

1.26 

1.37 

1.38 

1.17 

1.08 

t.01 

For  these  amplification  factors,  the  following  bounds  for  the  response 
spectrum  are  computed: 

A'  =  A  X  PGA  =  2.71  X  0.^  =  1.3!g 

V'  *  V  X  PGV  =  2.30  X  61  cm/sec  =  \A0cm/sec 

1/  -  D  X  PGD  =  2.01  X  45  cm  =  90cm 

The  resulting  elastic  response  spectrum  is  that  shown  in  Hgure  15.  As 
shown  there,  Newmarit  and  Hall  (1982)  connect  the  SDOF  acceleration  at 
8  cps  to  the  peak  ground  acceleration  at  33  cps,  thereby  completing  the 
high-frequency  portion  of  the  spectrum.  They  do  not,  however,  specify  at 
what  frequency  the  SDOF  displacement  should  be  connected  to  the  peik 
ground  displacement  to  complete  the  low-frequency  pcntion  of  the  spectrum. 
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In  itadies  and  analyses  separate  frmn  those  of  Newmark  and  Hall,  the 
Applied  TBchnology  Council  (ATC)  developed  another  form  of  the 
broiul-band,  smooth-shaped  response  spectra.  The  results  of  their  anal¬ 
yses  are  presented  in  publication  ATC  3-06.  dated  1978. 

The  ATC  3-06  broad-band,  smooth-shaped  spectra  were  based  on  the 
smoothed  mean  spectral  shapes  from  the  study  by  Seed,  Ugas,  and  Lysmer 
(1976)  shown  in  Figure  16.  These  average  spectra  are  based  on  104  re¬ 
cords  mostly  from  earthquakes  from  the  western  part  of  the  United  States 
and  having  a  Richter  magnitude  range  of  S.2S  to  7.5. 


Figure  16.  Average  acceleration  spectra  (Seed.  Ugas.  Lysmer  1976) 


Like  the  Newmark  and  Hall  analysis,  the  end  product  of  the  ATC  analy¬ 
sis  is  a  graph.  The  ATC  graph  is  a  linear  graph  on  which  the  periods  of 
SDOF  systems  are  plotted  against  the  normalized  spectral  acceleration 
(pseudoacceleration  divided  by  the  peak  ground  acceleration).  An  exam¬ 
ple  of  the  ATC  3-06  response  spectra,  drawn  for  5  percent  damping,  is 
shown  in  Figure  17;  a  distinct  advantage  of  the  ATC  response  spectra  is 
that  the  various  soil  conditions  are  included  directly  with  the  spectra. 

In  Figure  17,  the  nonsite-specific  response  spectra  appropriate  for  use 
at  a  given  site  are  distinguished  between  the  three  site  classifications 
given  in  Table  2. 

The  abscissa  of  the  ATC  3-06  response  spectra  of  Figure  17  is  normal¬ 
ized,  that  is,  it  is  the  ratio  of  peak  SDOF  pseudoacceleration  to  peak 
ground  acceleration.  Note  also  that  Figure  17  is  a  linear  plot  rather  than  a 
log  plot. 
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Figure  17. 

ATC  3<0€  response  spectra,  5  percent  damping 

Tabl«2 

SoHTypM 

Tlfpe 

OMOrifNlOfW 

1 

A  sol  prone  wHh  ettier 

a.  AracMlwmnerielctieraciertiedbyeeheerwa»e»eloctty 
orealor  ttwn  2.500  Ipe  or  by  oSier  euRable  means  ol  ctaesMceSon. 

or 

b.  SSir  or  dense  soRoondHionwbaraaie  soil  depSi  is  leee  Own  200  fL 

2 

A  sol  proSto  wHh  dsnse  or  stiff  sol  oondWons,  where  «w  soi  depSi  exceeds 

200  tt. 

3 

A  sd  preIN  70  ft  or  mere  in  dspSi  snd  centaMne  mere  than  20 11  of  soft  to 
medhim-sM  day  but  not  more  than  40  ft  oT  soft  dsy. 

For  SDOF  systems  founded  on  rock  or  stiff  soils  (lype  1  soils),  die 
maxinmm  noniMlized  value  of  2  J  occurs  in  SDOF  systems  luiving  a  natural 
period  between  0.2  to  0.4  sec.  For  SDOF  systems  founded  on  lype  I  soils 
but  which  have  periods  longer  than  0.4  sec,  the  normalized  values 
decrease  in  amplitude  inversely  proportional  to  the  period  T,  as  indicated 
in  Rgnie  17  (recall  that  T  «  I/f)- 
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For  toiti  loflw  Ham  Type  1  soils,  the  peak  nonnalized  acce^irtioB  is 
ekMfiaed  tignifieaBtly,  exteediBg  to  periods  of  0.6  sec  for  Type  2  soils 
and  to  0.9  sec  for  Type  3  soils.  Even  when  acceleration  does  start  to  de¬ 
crease  in  Types  2  and  3  soils,  the  rate  of  decrease  is  less  than  that  in  stif- 
fer  soils.  Such  long  periods  are  more  characteristic  of  flexible  structural 
systems  than  rigid  structural  systems,  indicating  that  when  soils  are  soft, 
spectral  accelerations  are  likely  to  be  high  for  both  flexible  and  rigid  struc 
tural  systems. 

The  greatest  dynamic  amplification  of  ground  motion  occurs  in  struc¬ 
tures  when  the  fundamental  period  of  the  structure  is  close  to  the  charac¬ 
teristic  period  of  the  ground  motions.  This  means  that  flexible  or 
long-period  structures  on  soft  sites  will  respond  more  to  earthquake 
ground  motions  than  flexible  structures  on  stiff  sites.  The  ATC  3-06  stan¬ 
dard  spectra  for  various  site  condidoos  attempts  to  capture  this  phenomenon. 

As  a  matter  of  interest,  the  smoothed  plot  for  rock  foundations  as  it 
was  proposed  in  the  Seed.  Ugas.  Lysmer  study  is  shown  in  Figure  18.  It 
has  already  been  noted  that  the  ATC  curves  evolved  from  this  study. 
Shown  also  for  comparison  are  the  plot  as  it  was  finally  adopted  for  the 
ATC  3-06  curve  and  the  corresponding  plot  taken  from  a  Newmark  and 
Hall  response  spectra  for  the  same  rock  foundation. 


Figure  18.  Comparison  of  design  response  spectra  for  rock  foundations 


Comparing  the  ATC  curve  with  the  Seed  curve  (a  “mean”  curve)  shows 
that  ATC  truncated  Seed’s  peak  normalized  acceleration  at  2.5.  This  maxi¬ 
mum  of  2.5  is  qiplied  to  all  soils  in  the  ATC  spectra  (see  Figure  17),  not 
just  rock.  Also,  ATC  more  than  tripled  the  range  of  periods  which  would 
experience  the  maximum  normalized  acceleration  of  2.5. 

Comparing  the  two  major  response  spectra  of  Newmark  and  Hail  and 
ATC  3-06  shows  that  the  maximum  magnitude  of  the  ATC  response  is 
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significantly  higher  than  that  of  Newtnark  and  Hall  (a  “median”  curve), 
and  remains  higher  for  all  values  of  T.  It  is  concluded  that  the  ATC  3-06 
spectrum  is  somewhat  more  conservative  than  the  Newmark  and  Hall  spec¬ 
trum  for  all  frequencies. 

TIm  ATC  3-06  spectra  also  appear  in  the  1991  Uniform  Building  Code 
as  a  general  nonsite-specific  design  requirement.  They  are  included  as 
well  in  the  recommended  Lateral  Force  Requirements  of  the  Structural 
Engineers  Association  of  California  (SEAOC)  and  are  also  used  in  the 
tri-services  manual  TM  5-809-10-1.  SEISMIC  DESIGN  MANUAL  FOR 
ESSENTIAL  BUILDINGS.  The  ATC  3-06  IVpe  I  response  spectra  are 
used  in  EM  1 1 10-2-2401  intake  tower  example  problems. 

Standard  spectra,  such  as  ATC  3-06.  are  intended  to  be  used  only  for 
preliminary  structural  evaluations  when  site-specific  response  spectra  are 
not  yet  available  for  the  site.  Another  procedure  for  developing  standard 
spectra  to  be  used  in  preliminary  structural  evaluations  is  described  in 
ER  1110-2-1806. 
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3  Distributed-Mass  Systems 


Period  of  Oscillation  of  Distributed-Mass 
Systems 

A  tnie  SDOP  system  is  one  in  which  only  one  deflection  pattern  can 
occur.  All  of  the  Inmped-mass  systems  shown  in  Figure  1  are  such  SDOF 
systems.  In  solutions  to  SDOF  systems  there  was  never  any  need  to  con¬ 
sider  multiple  patterns  of  oscillations. 

A  distributed-mass  system  such  as  the  tower  shown  in  Figure  19  is  not 
an  SDOF  system.  The  tower  could  have  a  number  of  patterns  of  oscilla¬ 
tions  as  shown,  and  it  could  also  oscillate  in  two  or  more  of  these  patterns 
at  the  same  time.  The  distributed-mass  system  of  Figure  19  is,  in  fact,  a 
multiple-degree-of-freedom  (MDOF)  system  with  an  infinite  number  of 
degrees  of  freedom. 


A4 


If.  however,  it  is  assttmed  that  the  tower  oscillates  in  tmly  a  single  de¬ 
flection  shape  at  any  one  time,  the  system  could  be  analyzed  as  an  SDOP 
system.  Such  an  assumption  is  the  basis  of  analysis  of  distributed-mass 
systems,  i.e..  that  the  system  can  in  fact  be  limit^  to  a  single  pattern  of 
deflections  at  a  particular  time.  Hie  Rayleigh  method  provides  such  an 
analysis  of  distributed-mass  systems. 

In  the  Rayleigh  method,  it  is  assumed  that  the  deflection  u(x,t)  is  lim¬ 
ited  to  a  single  shape  and  that  the  single  shape  can  be  expressed  as  vari¬ 
ables  separable: 

uOcj)  =  yfixMt)  (38) 

Inherently,  the  variables-separable  solution  uncouples  the  rulution  for 
displacements  from  the  solution  for  frequency:  each  can  be  solved  inde¬ 
pendently  of  the  other.  The  time  Junction  v(t)  defines  only  the  magnitude 
of  the  horizontal  displacement  at  some  reference  level  r  at  any  time  t.  The 
displacement  function  y  becomes  only  a  shape  function,  defining  the 
relative  shape  of  the  tower  at  any  point  x  along  its  height  as  shown  in 
Figure  20. 

For  the  sake  of  simplicity,  the  reference  level  r  in  this  report  is  always 
taken  at  the  top  of  the  cantilever,  that  is.  x  in  Equation  38  equals  L. 


Figure  20.  Shape  function 


Chapters  Dtotrtbuted-Mam Systems 


OMCBT  KB  of  tihe  SDOT  tytttm  of  Pignie  20  is  given  by: 
KBm  Mt2 

«  V4  m(r)  [uMl^dx 
0 

The  stored  potential  energy  PE  in  the  system  is  the  strain  energy  due  to 
flexiue: 


\A  r 


where  M  is  the  moment  mi  the  cross  section  and  6  is  the  rotation  of  the 
cross  section  about  the  neutral  axis.  The  rotatimi  6  is  given  by  the  Ber- 
nonlli  equation; 

a  L  M.  ih. 

p  “  H  ’ 


The  final  form  of  the  PE  equation  is  then: 

P£  =  f  £f(*) 

0 


(42 


An  important  observation  is  appropriate  at  this  point.  Note  that  the 
potential  energy  used  in  the  Rayleigh  solution  is  limited  to  the  flexural 
moment-times-iotation  strain  energy.  Shear  deformations,  axial  deforma¬ 
tions,  and  torsional  rotations  are  all  ignored.  The  Rayleigh  solution  is, 
therefore,  most  accurate  in  those  applications  in  which  the  flexural  energy 
dominates  all  other  strain  energy  sources  to  such  an  extent  that  they  may 
be  safely  ignored. 

It  should  be  noted  that  for  squat  towers  the  shear  displacement  can  be  a 
signiticant  part  of  the  total  lateral  displacement  (flexural  and  shear).  In¬ 
take  towers  with  height-to-width  aspect  ratios  of  3  or  less  should  be  de¬ 
signed  for  earthquake  ground  motions  using  procedures  that  include  shear 
stiffness  cqwbility. 

The  variation  in  KE  at  any  time  interval  must  be  equal  and  opposite  to 
the  variation  in  the  PE  over  the  same  interval  since  there  is  no  damping  in 
this  conservative  system.  The  resulting  solutions  of  Equations  39  and  42 
are  rather  complicated  solttions  in  the  calculus  of  variations  (Paz  1991} 
which  yields  the  Rayleigh  solution: 

m*ii(0  +  **v(0  =  (43 
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m*  «  J  iiiix)y/^dx  •  gtneraUzed mass 

0 

(44a) 

4;*  «  J  £J(x)  [y^dx  »  generalized  stiffness 

0 

(44b) 

-iig  m(x)  qafx  «  generalized  load 

(44c) 

In  Equatioos  43  and  44,  the  asterisks  denote  that  the  functions  are  general- 
iud  fiinctiaos  which  occur  at  the  pefetOKe  fevel  r  at  die  top  of  the  cantilever. 

A  carefhl  examination  of  Equatims  43  and  44  reveals  that  the  Rayleigh 
solution  is  actually  a  solution  for  a  lumped*ma8s  system  rather  than  fmr  a 
distributed-mass  system.  Such  an  equivalent  system  is  shown  in  Fig¬ 
ure  21. 

As  indicated  in  Hguie  21.  the  Rayleigh  solution  given  by  Equation  43 
is  an  equivalent  solntimi  for  the  oscillations  of  a  generalised  mass  m*  at 
die  reference  level  r  with  respect  only  to  time.  The  solutimi  for  i>(t)  is  no 
longer  a  function  of  the  height  x.  This  indicates  that  the  Rayleigh  solutirm 
transfmms  an  MDOF  distrilmted-mass  system  into  an  equivalent  SDOF 
lunqied-mass  system  at  the  reference  level  r. 


Figure  21 .  Equivalent  Rayleigh  iumped-maas  system 
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BfMikm  43  dost  aot  bwtode  dtaipUif.  i«.,  theie  u  oo  loM  of  eaeify 
nyiriMO  ia  tto  qrateoB.  It  was  Bhowa  e«rUer,  however,  that  the  natnnl 
period  of  oacillatioo  is  aot  affected  by  <fauiq>iag  (for  low  fevels  of  danq;^ 
ia|),  so  the  period  cooipated  by  ase  of  Bqaation  43  will  still  be  correct 

The  sointioa  Equatioa  43  takes  the  form: 

ikt)  +  »  -(V"**)  (^5) 

where 


it  (46) 

"  0 

The  natural  frequency  of  this  system  remains,  of  course,  the  unforced  fre¬ 
quency.  Again,  the  unforced  solutimi  for  i>(t)  is  harmonic.  i.e.. 

vit)  s  Dq  sin  cor  (47) 

and 


iM  a  >|^(*)Oo  sin  <nr  (48) 

where  Vq  is  the  maximum  value  of  displacement  at  the  reference  level  r. 

With  no  danq>ing.  the  energy  in  this  conservative  system  remains  con¬ 
stant.  The  sum  of  K£  (Equation  39)  and  PE  (Equation  42)  is  therefore  a 
constant.  At  maximum  amplitude,  the  velocity  is  zero,  and  all  the  energy 
in  the  system  at  that  time  is  in  the  form  of  PE: 


PE 


•a 


£/(*)  t^dx 


(42) 


Equation  48  is  substituted,  yielding: 
PE  =  ‘AoJ  £i(X) 


(49a) 


or 


PE=Viy^k*  (49b) 

where  Vq  is  the  displacement  at  the  reference  level  r. 

Similarly,  when  displacement  u(x,t)  is  zero,  all  of  the  energy  in  the  sys¬ 
tem  is  in  the  form  of  KE: 
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Jffi-  Vi  f  m(;c)  luOcj)]^dx 


(39) 


Equation  48  is  substituted,  where,  with  cos((Dt)  =  1  when  i(<x,t)  is  at  its 
maxinnini: 


KE 


m(jc)(v(*)l^dx 


(50a) 


or 


Vit^oAii*  (50b) 

Equations  49  and  50  must  be  equal  in  a  conservative  system.  Equating 
them  produces  the  result: 

a  k*/m*  (5ia) 

or  equivalently, 

r  =  2*  Vm*/3k*  (51b) 

When  Equation  51a  is  substituted  into  Equation  45,  the  result  is  the 
final  form  of  the  Rayleigh  solution: 

ikO  +  wMO  » -(V"**)V^> 

where  the  natural  period  T  and  the  circular  frequency  oo  are  those  given  by 
Equations  51a  and  b.  _ 

The  factor  Lg/m*  on  the  right-hand  side  of  Equation  52  is  an  important 
feature  of  the  Rayleigh  solution.  It  is  a  dimensionless  number,  a  multi¬ 
plier  that  is  a{q>lied  to  the  ground  acceleration  ag(t).  As  such,  it  becomes 
a  part  of  the  forcing  function  and,  consequently,  it  automatically  becomes 
the  same  multiplier  of  the  spectral  quantities  Si>,  Sy,  and  (to  be  dis¬ 
cussed  later). 

Equation  52  shows  that  the  factor  Lg  by  itself  is  the  "normalization’* 
factor  between  the  actual  ground  acceleration  Og(t)  and  the  response  v  by 
the  generalized  mass  m*.  When  divided  by  m*,  the  entire  factor  Lg/m*  be¬ 
comes  the  normalization  ratio  factor  between  Qg(t)  and  i>  per  unit  of  gener¬ 
alized  mass. 

Equation  51a  has  more  than  one  solution.  The  only  requirement  for  a 
solution  is  that  die  strain  energy  stored  in  a  given  deformation  shape  be  ex- 
acUy  equal  to  the  work  done  by  the  mass  in  moving  into  that  shape.  For  a 
distributed  mass,  there  are  an  inHnite  number  of  shape  functions  that  will 
satisfy  all  the  energy  relations  and  thus  be  a  solution  to  Equation  51a. 
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Eaeli  tncb  thi^ie  ftmction  will  have  its  own  corresponding  circular  fre¬ 
quency  os  and  period  T. 


In  order  to  compute  the  period  of  the  structure  from  Equation  51b,  it  is 
necessary  to  know  the  mass  m(x)  and  its  distribution,  the  stiffness  EI(x) 
and  its  distribution,  and  the  shape  function  y.  Where  El  is  constant  or  var¬ 
ies  in  discrete  steps,  it  is  known  from  the  Bernoulli  equation  that  the  re¬ 
storing  force  on  the  system  is  the  fourth  derivative  of  the  displacement 
function  y: 

^  £_(m\  ^  (53) 

This  restoring  force  at  any  point  is  equal  to  the  spring  constant  k  times 
the  deflection  y,  with  the  direction  of  the  force  always  opposite  to  the  di¬ 
rection  of  the  ^flection  y: 


*¥  =  -P  (54) 

Equating  these  two  values  for  p  yields  an  equation  of  the  form: 

y""  =  (constant)  (y)  (55a) 

or, 

y""  -  (constant)  (y)  =  0  (55b) 

The  solution  of  an  equation  of  this  form  is  given  by: 

y  =  A  sin  (az)  +  B  cos  (az)  +  C  sinh  (az)  +  D  cosh  (az)  (56) 

where  z  =  x/L  and  a  is  a  constant  involving  the  rotational  frequency 
(0  (a*  =  mjs^/El,  Clough  and  Penzien  1993). 

There  are  only  four  constants  of  integration  in  the  solution,  which  lim¬ 
its  the  boundary  conditions  to  four.  There  are  frve  unknowns,  however, 
since  a  is  also  unknown.  Therefore,  a  fifth  condition  is  required.  One 
way  to  develop  the  fifth  condition  is  to  assign  some  arbitrary  value  (such 
as  unity)  to  y  at  the  reference  level  r  (at  the  top  of  the  cantilever,  see  Fig¬ 
ure  21).  The  solution  will  then  be  a  reference  solution,  or  normalized  solu¬ 
tion,  in  terms  of  unit  displacement  at  r. 

An  example  will  illustrate  the  format  of  the  solution.  The  first  three 
shape  ftinctions  for  the  cantilever  of  Figure  22  will  be  found.  The  four 
boundary  conditions  for  the  cantilever  are: 

1)  The  deflection  at  the  base  of  the  cantilever  is  zero 
atzaO,  y  =  0 
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jure  22.  Typical  distributed-mass  cantilever 

2)  The  slope  of  the  elastic  carve  at  the  base  is  zero 
atzvO,  ifr'aO 

3)  The  moment  at  the  free  end  is  zero 
at  z  =  I,  *  0 

4)  The  shear  at  the  free  end  is  zero 
at  z  «  1,  ly"' «  0 

The  derivative  forms  of  Equation  S6  are: 
qr  «  A  sin  (az)  +  B  cos  (az)  C  sinh  (az)  -i-  D  cosh  (az) 

V  °  a  A  cos  (az)  -  aB  sin  (az)  +  aC  cosh  (az)  +  aD  sinh  (az) 
yf''  K  -a^A  sin  (az)  -  a^B  cos  (az)  +  a^  sinh  (az)  +  a^D  cosh  (az) 
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fr  .t^;i«^^;^3MpfHl4t%iii(«()4>a%«Mli(ai>'i>a’Dttaii(as) 

»  ft^Aiil  (tt)  *  »*B  cot  (as)  ♦  t^  ciali  (as)  -¥  a^  coah  (as) 

Tkeae  an  tadMftniad  iirto  the  boandaiy  cimditioas  to  find: 

Ats«Oki|r  «0*B4>D 

Ats^O.if'  »0»  A-fC 

At  z  ■  ],  Iff''  at  0  a  -A  tin  (a)  >  B  cos  (a)  -i-  C  sinh  (a)  -i-  D  cosh  (a) 

At  z  a  1,  Y**  >■  0  a  -A  cos  (a)  B  sin  (a)  -f  C  cosh  (a)  +  D  sinh  (a) 

A  fifth  conditioo  is  now  imposed  to  <rfitain  the  reqaiied  fifth  equation. 
The  vahM  qr  at  the  free  end  of  the  cantilever  (reference  level  r)  is  as- 
tified  a  vatne  of  unity,  yielding: 

At  z  a  1,  tp  a  i  M  A  sin  (a)  B  cos  (a)  *  C  sinh  (a)  D  cosh  (a) 

A  set  of  equations  such  as  these  five,  havmg  a  variable  coefficient  in 
the  matrix,  is  called  a  characteristic  problem  or  an  eigenproblem.  A  solu¬ 
tion  win  exist  only  for  particular  values  of  the  variable  a;  these  values  of 
a  are  called  the  characteristic  values  or  the  eigenvalues. 

There  are  an  infinite  nundnsr  ct  values  of  a  in  any  distributed-mass  sys¬ 
tem.  Each  value  of  a  corresponds  to  a  functitm;  the  lowest  value  of 
a  corraqKmdiag  to  the  first  shape  function,  the  next  lowest  value  of  a  cor¬ 
responding  to  the  second  shqie  function,  and  so  on.  For  the  sake  of  brev¬ 
ity,  only  three  shiqie  functions  are  presented  here. 

These  ft ve  simultaneous  equaticms  are  solved  using  conventional  com¬ 
puter  strflware  (e.g..  Maple,  MathCAD,  Matlab,  Matbematica,  etc.), 
diough  there  are  also  specialized  hand-held  calculators  available  (e.g., 

HP  48S,  TI 81,  n  85,  etc.)  that  could  solve  such  a  set  of  equations.  The 
first  three  solutions  are: 


a-  i.araio 

a-  4.60409 

a-  7.85476 

A-40.36708 

Axasoozs 

A--»0.4e961 

1 

■ 

a 

B--»OJOOOO 

B--O.S0000 

! 

■ 

o 

C--»OJOe23 

C- -0.49961 

D-4OJ0000 

0- -0.60000 

D-40.90000 

Dqpmiding  upon  the  particular  software,  the  trial  values  (rf  a  are  entered 
ittimuctively  by  the  engineer.  The  lowest  value  of  a  for  which  a  solution 
am  be  obtdned  will  yield  the  ftrst  sluqm  function.  The  next  lowest  value 
of  a  will  yield  the  seomd  shape  function,  and  so  on.  It  is  recommended 


•ad  so  on.  It  is  recominended  that  the  gr^hics  cafMdiility  of  the  st^ware 
be  nsed  to  verify  each  shape  fonction  so  obtained.  The  number  of  curva¬ 
tures  correqwnds  to  the  number  of  the  shape  function  (see  Figure  23). 

Therefore,  the  equations  of  the  first  three  shape  functions  are: 

-  0.367  sin  (1.87Sz)  -  0.S00cos  (1.87Sr) 

-  0.367  sinh  (1.875z)  -t-  0.S00  cosh  (1.87Sz) 

Vj  *  -0-^  sin  (4.694Z)  +  0  JOO  cos  (4.694z) 

0.309  sinh  (4.694z)  -  0.500  cosh  (4.694z) 

^3  «  0.500  sin  (7.855z)  -  0.500  cos  (7.855z) 

-  0.500  sinh  (7.855z)  -t-  0.500  cosh  (7.855z) 

The  configurations  of  the  shiqw  functions  yf  corresponding  to  these  first 
three  solutions  are  shown  in  Rgure  23. 


Figure  23.  Shape  functions  lor  a  cantiiever  having  a  uniformiy  distributed 
mass 


The  three  equations  for  y  are  now  used  to  evaluate  the  generalized 
functions  for  m*.  k*.  and  Lg  given  by  the  integrals  of  Equations  44a  and 
b  and  Equation  46  (or,  equivalently.  44c).  The  integrals  may  be  evaluated 
by  conventional  software  or  by  specialized  hand-held  calculators.  The 
generalized  values  for  the  three  shape  functions  of  Figure  20  are  listed  in 
Table  3,  along  with  the  rotational  frequencies  and  the  natural  periods. 
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TiMS 

OtmrallMdVtfii— foraPtotflbutedllBfCwtllwr 

Ihaps  FuwsMsii  l 

•kape  nmeUoN  2 

Skepe  FuneUen  3 

m*  >038000  eg. 

m*- 0.24800  ng. 

m*  -  008018  ng. 

k*  .  3.08088  EM.* 

■(*-121.3083  EM.* 

■(*>8820417  EM.* 

L,  -  0.38148  mjL 

L,  .  -0.21701  ng. 

■^- 0.12758  ng. 

La/m*- 1.86888 

Lgm*  - -006811 

igm*  -  0.50885 

«/>  12.30  (EMig.*) 

tt*  -  488.5  (Eimg.*) 

m  3808  {BVmjL*) 

Tminfm 

T-2k/« 

T  m  7xtm 

Flexural  Stiffness  El  for  OBE  and  MDE 

The  value  of  k*  given  in  Table  3  is  a  Ainction  of  the  modulus  of  elastic¬ 
ity  B  and  the  moment  of  inertia  I.  For  concrete  structures  under  the  OBE, 
the  value  of  El  is  taken  conservatively  as  that  of  the  uncracked  concrete 
section.  The  uncracked  value  of  El  is  used  for  determination  of  both  the 
accelerations  and  displacements  under  an  OBE  earthquake. 

For  the  maximum  design  earthquake,  it  is  recognized  that  the  worst- 
case  values  of  acceleration  and  velocity  will  occur  if  the  section  does  not 
crack  this  because  the  higher  the  stiffness,  the  higher  the  structural  system 
frequency  and  corresponding  spectral  acceleration  and  velocity.  The  un¬ 
cracked  section  is,  therefore,  used  when  calculating  accelerations  and 
velocities  under  the  MDE.  For  displacements,  however,  the  largest  dis¬ 
placements  will  occur  when  the  concrete  cracks,  resulting  in  a  reduction 
in  stiffiiess  which  increases  the  period  of  the  structure  and  its  correspond¬ 
ing  spectral  displacement.  Under  the  MDE,  therefore,  the  uncracked 
value  of  El  is  used  to  determine  accelerations.  To  determine  displace¬ 
ments,  a  flexural  stiffness  equal  to  half  the  gross  section  flexural  stiffness 
is  generally  used.  This  value  is  considered  an  approximate  typical  effec¬ 
tive  stiffness  for  reinforced  concrete  structures  loaded  to  yield  levels. 


Combined  Effects  In  DIstrlbuted-Mass  Systems 

The  values  of  the  period  T  found  from  the  Rayleigh  solution  are  used 
to  ffnd  the  spectral  values  of  Sd,  Sy,  and  S^,  either  from  standard  spectra 
such  as  Newnnaric  and  Hall  or  ATC  3-06  or  from  a  design  response  spectra 
developed  for  the  specific  site.  There  will,  of  course,  be  three  sets  of  spec¬ 
tral  values,  one  for  each  shape  function.  Since  the  cantilever  can  experi¬ 
ence  all  three  (or  even  more)  of  the  oscillation  shapes  at  the  same  time. 
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the  spectral  values  must  now  be  combined.  A  means  to  combine  the  spec¬ 
tral  values  from  several  shape  functions  is  presented  in  this  section. 

Finding  a  rational  means  to  combine  several  shape  functions  is  a  major 
consideration  if  a  set  of  SDOF  solutions  is  to  be  taken  as  the  solution  to 
an  MDOF  problem.  Because  a  system  may  oscillate  in  more  than  one  fre¬ 
quency  at  any  one  time,  a  combination  of  SDOF  shapes  will  undoubtedly 
occur  in  any  MDOF  oscillating  system.  For  each  of  these  participating 
shapes,  a  spectral  maximum  displacement  Sq  can  be  found,  along  with  its 
corresponding  spectral  maximum  pseudovelocity  Sy  end  its  spectral  maxi¬ 
mum  pseudoacceleration  The  spectral  values  are  the  maximum  values 
for  each  frequency.  Since  it  is  unlilmly  that  all  these  maxima  will  occur 
at  the  same  time,  some  reasonable  basis  for  combining  these  maximum 
values  is  obviously  needed. 

There  are  two  methods  in  common  use  for  combining  several  SDOF  so¬ 
lutions  to  obtain  a  single  MDOF  solution: 

•  square  root  of  the  sum  of  the  squares  (SRSS),  and 

•  complete  quadratic  combination  (CQC) 

The  older  and  more  universal  of  these  two  methods  is  the  SRSS  method, 
which  is  an  approximate  method.  Any  number  of  spectral  values  can  be 
combined  by  SRSS,  but  the  usual  number  is  three.  When  combined  using 
SRSS,  the  spectral  values  for  n  SDOF  solutions  are: 


So  =  ^  <57a) 

Sy  =  V4,  +  4  +  +. . Si,  (57b) 


The  more  recent  CQC  method  is  based  on  random  vibration  theory  and 
it  accounts  for  the  interaction  between  two  shape  functions  that  have  peri¬ 
ods  relatively  close  to  each  other.  Hie  method  uses  a  cross-function  coef- 
frcient  Pij,  where  i  and  j  are  any  two  shape  functions  and  where  r  =  Tj/Tj: 

=  (1  +  r)  (58) 

”  (1  -  +  4pV  (1  +  r)^ 

Figure  24  shows  a  graph  of  py.  For  values  of  damping  less  than  5  per¬ 
cent,  the  periods  must  be  within  about  25  percent  of  each  other  for  signifi¬ 
cant  cross-function  interaction  to  occur. 
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Figura  24.  CroM’functkm  coefficient  py 


The  spectral  values  for  a  CQC  combination  are  given  by: 

VN  H 

X  X 

i-iy-i 

VTi — n - 

X  X  ^vr 

p — R - 

-  ■'•X  2 

<-U-I 

When  the  periods  between  functions  are  well  separated,  the  cross-function 
coeffmient  q)proaches  zero  for  all  cases  where  i  w  j.  For  the  case  i.’ 
which  i«j,  Ae  value  of  Pj.  becomes  1,  and  the  CQC  method  degenerates 
into  the  SRSS  method.  Most  large  structures  will  have  distinct  separation 
between  periods  and  will  fall  into  this  default  category  in  which  tte  SRSS 
method  will  apply. 

Table  3  shows  that  the  normalization  factor  Lg  can  have  different  alge¬ 
braic  signs  for  different  siuqm  functions.  When  the  CQC  method  is  used, 
diose  signs  must  be  preserved  throughout  the  calculations.  Further,  when 
combining  east-west,  north-south  and  up-down  contributions  to  any  spec¬ 
tral  quairtity,  the  contributions  must  be  added  algebraically  before  being 
used  either  in  the  SRSS  method  or  the  CQC  method  (EC  1 1 10-2-60S0). 

The  following  example  illustrates  the  procedure  for  finding  the  spectral 
response  of  a  distributed-mass  system  in  which  Rayleigh’s  solution  is 
used  to  find  the  period  T. 


(59a) 

(59b) 

(59c) 
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Example  Solution  for  a  Stopped  Tower 


ProbiMi  •tttMiMnt 

Calculate  the  forces,  overturning  moment,  and  displacements  on  the 
square  step-tapered  intake  tower  under  the  maximum  design  earthquake, 
using  the  uncracked  section  for  computing  forces  and  the  fully  cracked 
and  rotated  section  for  computing  die  displacements.  The  Type  1  ATC  3-06 
response  spectrum  is  to  be  used.  Maximum  ground  acceleration  is  0.45g, 
the  foundation  is  rock,  Y^onc  ^  Ib/ft^.  and  Gross  =  3.6  x  10^  Ib/in.^ 


CoiiHHitt 

a.  Natural  periods  of  oscillmion  for  the  first  three  shape  functions. 

b.  Spectral  displacements,  velocities,  and  accelerations  for  the  first 
three  shape  functions. 

c.  Combined  acceleraticms  for  the  first  three  shape  functions. 

d.  Base  shear  under  the  combined  accelerations. 

e.  Overturning  moment  under  the  combined  accelerations. 

/  Displacements  of  the  cracked  and  damaged  tower. 


-  1984  slugs^ft  =  l.OOm^ 

/3  =  b}?/\2  =  72,470yr*  =  1.000/3 

=  2534  slugs^ft  -  1.28m^ 

I2  =  bl^/\2  =  106,260^  =  1.466/3 

m^j  =  3261  sluga^ft  -  1.64m^ 

/j  =  b}?/\2  =  145.830yjf*  =  2.012/3 
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Btww  th«  ihapt  ftmctkw  is  diseoiHiMOQS  at  the  steps,  it  is  defined 
by  thrae  dbeoatiMKMS  eqaadnas; 

0  jfc  z  ^  V6:  yfm  AjSin  (az)  +  Bjcos  (az)  'f  CjSinh  (az)  +  Djcosh  (az) 
V6  j£  z  1  yfm  AjSin  (az)  -f  Bjcos  (az)  -t-  C2sinh  (az)  +  D2Cosh  (az) 
%  S  z  22  1:  ly  »  AjSin  (az)  *  Bjcos  (az)  +  C3sinh  (az)  *  DjCOsh  (az) 
Boundary  conditions  are: 

At  z  «  0,  Yj  s  0  and  Vj  «  0 
At  z  ■  Yi  ■  ¥2;  ¥^1  *  Vi*'  =  M2  and  Vj  =  V2 
At  z  «  H  ¥2  ■  ¥3;  V2  “  ¥'3*  ^^2  “  *^3  “<*  ^2  *  ^3 
Atzvl,  MjsOandVjsO 
when 


M  ■  immient 


V  sshear 

These  12  boundary  conditions  apply  to  13  unknowns.  The  IS***  condi> 
tion  is  obtained  by  setting  the  displacement  Y3  »  1  at  the  top  of  the  tower. 

At  z  «  1.000.  Y3  "  i  ooo 

The  three  equations  of  the  shape  function  an  substituted  into  these 
13  boundary  conditions,  producing  13  equations  in  13  unknowns.  The 
equmions  an  independent  of  mass. 


0  at  B|  -f  D| 
Om  A|  -f  C| 


0  »  (AJ-A2)  sin  (V6a)  +  (B.-B-)  cos  (^a) 

+  (<-j-C2)  (Vja)  +  (D1-D2)  cosh  (Vta) 

0  •  (Aj'Aj)  cos  (VSa)  -  (B1-B2)  sin  (!6a) 

+  (C,-C2)  cosh  (V^)  +  (D,-D2)  sinh  (V6a) 


0  «  -(I|A,>l2A2)  sin  (Vto)  +  (LBi-LB.)  cos  (l6a) 


T  *2r2f 

+  (I,C,-l2C2)  sinTi  G6a)  +  (liD,-l2D2)  cosh  (V6a) 


0  ■  -(IjAj-LA,)  co»  (V6a)  -  (LB, -LB-)  sin  (Via) 

+  (IlCj-IjCj)  cosh  (i^)  +  (I,D,-l2D2)  sinh  (Via) 

0  -  (A--A3)  sin  (^6a)  +  (B-^-Bj)  cos  (^ia) 

+  (Cj-Cj)  sinh  (^)  +  (Dj-Dj)  cosh  (Via) 

0  «  (Aj-Aj)  cos  (Via)  -  (Bj-Bj)  sin  (Via) 

(Cj-Cj)  cosh  (%a)  +  (Dj-Dj)  sinh  (Via) 

0  ■  -(I2A2-I3A3)  sin  (Via)  +  (LB2-I3B3)  cos  (Via) 

+  (I2C2-I3C3)  sinh  (Via)  +  (I2D2-I3D3)  cosh  (Via) 

0  »  -(I2A2-I3A3)  cos  (Via)  -  (I2B2-I3B3)  sin  (Via) 

+  (I2C2-I3C3)  cosh  (Via)  +  (I2D2-I3D3)  sinh  (Via) 

0  a  -A3  sin  (a)  -  B3  cos  (a)  +  C3  sinh  (a)  +  D3  cosh  (a) 

0  s  -A3  cos  (a)  +  B3  sin  (a)  +  C3  cosh  (a)  +  D3  sinh  (a) 

1  s>  A3  sin  (a)  +  B3  cos  (a)  -t-  C3  sinh  (a)  D3  cosh  (a) 

It  is  well  to  observe  at  this  point  that  the  boundary  conditions  will  al¬ 
ways  include  the  four  end  conations  (at  the  base  and  at  the  top)  plus  four 
conditions  for  each  discontinuity  in  moment  of  inertia.  There  will  be  as 
many  of  these  interior  discontinuities  as  there  are  changes  in  El.  Each  dis¬ 
continuity  will  simply  add  four  more  equations  to  the  set  of  simultaneous 
equations.  Since  the  conqiuter  will  solve  forty  equations  as  readily  as 
four,  the  additional  equations  cause  no  more  labor  than  making  some  addi¬ 
tional  entries  in  the  mathematics  software. 

The  foregoing  equations  also  show  that  the  shape  function  of  any 
distributed-mass  MDOF  system  is  independent  of  the  distribution  of  the 
mass.  Additional  masses  may  be  added  or  deleted  at  will,  but  as  long  as 
the  stiffness  is  unchanged,  the  shape  of  the  deflection  curve  will  remain 
the  same.  The  mass  and  its  distritmtion  will  affect  the  magnitude  of  the 
displacements  (and  the  period  T),  but  not  the  shape  function. 

The  13  equaticms  of  the  example  are  solved  by  the  mathematics  soft¬ 
ware,  as  usual,  to  End  the  coefficients  at  the  three  steps.  The  coefficients 
are  listed  in  Table  4.  Trial  entries  of  2, 5,  and  8  for  the  value  of  a  were 
used  in  finding  the  three  shape  functions  shown  in  Table  4. 
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8lwp»  Function  CfMfllcImitt 

nraltlM|MF4iiiellon  SMond  SImim  FuncUen  TMrd  SIwp*  Function 

■  Bm4jnm  tmJMsn 

Flratstcp,0dzsi1/3 


The  equations  defining  the  first  shape  function  are: 

0  ^  2  ^  V6:  V  =  0.241  sin  (2.150z)  -  0.357  cos  (2.150z) 

-  0.241  sinh  (2.150z)  +  0.357  cosh  (2.150z) 

16  ^  z  ^  V  =  0.257  sin  (2.150z)  -  0.463  cos  (2.150z) 

•  0.400  sinh  (2.150z)  >•-  0.510  cosh  (2.150z) 

^  z  ^  1 :  y  =  0.235  sin  (2.150z)  -0.554  cos(2.150z) 

-  0.660  sinh  (2.150z)  +  0.757  cosh  (2.150z) 


CtwplirS  DteWbutid-iylaM  Syttom 


The  generalized  nuM  ni*  end  generalized  stiffness  k*  for  the  first 
function  ere  found  by  evelueting  the  integrals  of  Equetions  44e  end 
b.  For  the  uncracked  section: 

m»  m^L  (0.241  sin  (2.150z)  -  0.357  cos  (2.150z) 

~  0.241  sinh  (2.130z)  +  0.357  cosh  (2.150z)]^<fe 

+  mjjL  (0.257  sin  (2.150z)  -  0.463  cos  (2.150z) 

-  0.400  sinh  (2.150z)  +  0.510  cosh  (2.150z)l 


+ 


f  (0.235  sin  (2.150z)  -  0.554  cos  (2.150z) 
dsr 

-  0.660  sinh  (2.150^  +  0.757  cosh  (2.150z)]  ^dz 


s  90,560  s/ugs 


*♦  a  r  (El/L^  (-1.114  sin  (2.150z)  +  1.650  cos  (2.150z) 

0 

-  1.114  sinh  (2.150z)  +  1.650  cosh  (2.150z)]  *dz 

+  (-1.188  sin  (2.150z)  +  2.140  cos  (2.150z) 

-  1.849  sinh  (2.150z)  +  2.357  cosh  (2.150z)l  ^dz 

+  f*  (£/,/L^)  (-1.086  sin  (2.150z)  +  2.560  cos  (2.150z) 

-  3.051  sinh  (2.150z)  +  3.499  cosh  (2.150z)l2dz 

a  36.15  X  10^  Ib/fi 
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^  (2.1500  o  0JS7cot  (2.1500 

~  0.241  sinh  (2.1500  +  0.357  cosh  (2.15001  di 

+  [0.257  sin  (2.1500  -  0.463  cos  (2.1500 

-  0.400  sinh  (2.1500  +  O.SIO  cosh  (2.15001  di 

+  [0.235  sin  (2.1500  -  0.554  cos  (2.1500 

-  0.660  sinh  (2.1500  •«-  0.757  cosh  (2.1S0z)]  dz 

a  151.215  5/uss 

The  remaining  constants  are  calculated  from  the  following  values; 
Circular  frequency  m  =  "4k^/m*  -  19.98  rad/sec 

Period  T  =  ht/tn  =  0.314  sec 

Frequency  /  =  l/T  =  3.186  cps 

The  values  for  the  other  shape  functions  are  found  similarly. 


Flrat  Shape  FuneUon 

Second  Shape  Fiinetion 

Thiid  Shape  Function 

ffl*  m  90,S60  akiga 

m*  >  77,800  aluge 

in*«  79,700  tluos 

k*-36.2xl0^lbfft 

k*  -  851  X  10*  lim 

k*  -  6,116  X 10*  bm 

1^.151, 215  ahins 

L„- -79,480  ekigs 

L„  B  50,540  slugs 

IVm*- 1.670 

-1.022 

IVn*  -  0.634 

•  ■  19.96  raiMac 

«■  I05rad/Sec 

»m277  rad/Scc 

T- 0.314  see 

T-0.080eec 

T  B  0.023  sec 

The  spectral  accelerations  are  found  from  Figure  17  for  the  periods  just 
calculate. 


Fbet  Shape  PUncUen 

Second  Shape  Fimctlcn 

Third  Shape  FuneUcn 

Sa- 1.125  g 

8a  B  0.720  g 

Sa  b  0.554  g 

CtHpUrS  DtaWbutMl-MaM  Systwm 


Effects  of  Submergence 


It  is  again  noted  that  the  foregoing  solution  for  the  period  of  the  exam¬ 
ple  step-tapered  tower  is  mathematically  exact.  The  solution  in  common 
use  for  the  period  of  such  towers  when  they  are  used  as  intake  structures 
(EM  1110-2-2401,  or  Goyal  and  Chopra  1989)  includes  the  additional 
mass  of  water,  both  inside  and  outside  the  intake  tower,  that  is  being  accel¬ 
erated  by  the  motion  of  the  tower.  Figure  25  shows  a  schematic  of  such 
added  mass  for  two  pool  levels.  Even  though  the  Chopra  procedure  for 
finding  the  added  mass  of  water  is  approximate,  it  indicates  that  the  added 
mass  that  lies  below  the  waterline  can  sometimes  triple  the  effective  mass. 
The  overall  result  is  to  increase  markedly  the  natural  period  T  above  that 
of  a  “dry”  tower. 


Figure  25.  Added  mass  due  to  acceleration  of  water 


The  increase  in  mass  due  to  the  added  water  poses  no  particular  prob¬ 
lem  in  the  Rayleigh  solution.  The  added  mass  of  water  is  computed  in 
stepped  increments  as  prescribed  in  the  Chopra  procedure  and  is  then 
added  to  the  mass  of  the  tower.  It  is  not  necessary  that  the  steps  in  the 
water  mass  match  the  steps  in  the  tower  mass.  If  the  steps  do  not  happen 
to  coincide,  the  steps  in  water  mass  will  simply  create  some  new  dis¬ 
continuities  in  the  integration  for  m*.  As  noted  earlier,  however,  the 
shape  function  remains  unchanged  by  the  addition  of  hydraulic  mass. 

Since  the  added  mass  does  not  affect  stiffness,  the  generalized  stiffness 
k*  remains  constant  for  all  levels  of  water  both  inside  and  outside  the 
tower.  Since  the  period  is  proportional  to  the  square  root  of  the  mass  m* 
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er«a» 

■t*  for  wfcMM  of  wbawioaco  aad  thra  lecoleohttiag  T  fw  each  of 
tfMM  lovols.  H^no  M  it  a  plot  of  the  renht  oS  toch  cakolatioBs  for  the 
VKUofin  ttep-ti^eiod  tower. 


Figure  26.  Effects  of  submergenoe  on  natural  period 


The  magnitude  of  the  added  hydrodynamic  mass  is  a  fimction  of  the 
cross-sectional  geometry  of  the  tower  and  the  depth  of  submergence. 
Gradual  changes  in  the  tower  geometry  produce  gradual  changes  in  effec¬ 
tive  mass  and  stiflhess,  whereas  abrupt  changes  in  the  tower  geometry  pro¬ 
duce  abrupt  changes  in  mass  and  stiffness.  The  effect  of  taper  on  the 
generalize  stifAiess  is  quite  distinct.  The  periods  will  be  distinctly  differ¬ 
ent  (as  much  as  50  percent  different)  between  a  tower  that  has  a  uniform 
cross  section  and  one  that  has  a  tapered  cross  section,  even  when  the  mass 
is  the  same  (EM  1110-2-2401). 


Effects  of  Bridge  Structures 

Another  feature  of  intake  structures  that  has  not  been  included  in  the 
previous  discussions  is  the  existence  of  the  bridge  structure  that  is  com¬ 
monly  placed  at  the  top  of  intake  towers  to  permit  service  access.  Such  a 
bridge  structure  is  shown  in  Figure  27.  In  an  earthquake,  the  existence  of 
the  bridge  strocmre  can  change  the  effective  mass,  and  in  some  cases  can 
add  a  restraint  Aat  most  be  included  as  a  boundary  condition. 
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Figure  27.  Bridge  structure  on  an  intake  tower 


Unless  the  bridge  is  anchored  to  the  tower  by  pinned  or  fixed  connec¬ 
tions,  large  relative  displacements  between  the  bridge  and  tower  can 
occur.  These  displacements,  unless  accounted  for  in  the  design,  can  lead 
to  hammering  or  to  a  loss  of  the  bearing  support  at  the  bridge  seat.  The 
preferred  solution  (not  always  done)  is  to  tie  the  bridge  and  tower 
together,  forcing  the  two  elements  to  perform  as  a  single  system.  The  po¬ 
tential  for  damage  doe  to  hammering  or  to  loss  of  support  is  then  elimi¬ 
nated.  The.  type  of  bearing  supports  at  the  bridge  seats  will  determine 
whether  the  bridge  can  be  handled  simply  as  an  added  mass  or  whether  a 
separate  boundary  condition  will  be  necessary. 

In  many  intake  structures,  the  bridge  beam  bearings  are  pinned  in  one 
direction,  but  in  many  others  they  are  pinned  in  two  directions.  Further, 
the  bridge  beams  might  also  be  pinned  at  their  far  ends  at  the  supporting 
pier.  The  flexibility  of  the  tall  piers  is  enough  to  dissipate  any  thermal  ef¬ 
fects  on  the  doubly  pinned  beams.  With  restraints  such  as  these,  lateral 
forces  will  be  created  at  the  top  of  the  tower  due  to  the  bridge. 

For  earthquake  motions  transverse  to  the  longitudinal  axis  of  the 
bridge,  the  force  on  the  pinned  bearings  at  the  tower  is  the  inertia  force  of 
the  half  of  the  bridge  supported  by  the  tower.  The  additional  inertial  force 
may  be  included  in  the  Rayleigh  solution  by  converting  the  vertical  beam 
reaction  at  the  tower  to  an  equivalent  mass  distributed  over  the  height  of 
the  bridge.  The  step  in  mass  thus  created  adds  yet  one  more  discontinuity 
in  the  distribution  of  mass.  The  integration  for  the  generalized  mass  m*  is 
then  made  as  usual;  the  extra  discontinuity  simply  adds  another  term  to 
the  integration. 
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Wok  oaatimo  trauvwrM  to  the  axii  of  die  btldie,  the  eddhioii  of  the 
hfidge  mMt  hee  no  eflbet  on  the  die^ecemoit  Auwdoo  y  of  tower; 
edflheM  remaiiit  vncheogui.  Incofar  as  displacemrats  are  coMemed, 

<mly  the  magnitodes  of  the  displacemeots  are  changed  by  attaching  the  sep- 
amle  bridge  mass. 


For  eardiquake  motions  parallel  with  the  longitudinal  axis  of  the 
bridge,  the  force  on  the  pinned  bearings  at  the  tower  depends  on  the  type 
of  bearings  used  at  the  pier  at  the  far  end  of  the  beams.  If  the  far  bearings 
are  rollers,  the  force  created  at  the  tower  would  be  the  inertial  force  of  the 
entire  mass  of  the  bridge,  not  just  the  reaction  at  one  end.  The  recalcula¬ 
tion  m*  would  proceed  as  before  but  with  the  larger  added  mass.  If  the 
far  bearings  are  pinned,  however,  the  force  at  the  tower  would  be  the  iner¬ 
tial  force  of  the  entire  bridge  plus  the  force  due  to  the  spring  constant  k, 
of  the  adjacent  pier  as  it  is  deformed  laterally  by  the  motions  of  the  tower. 
The  existence  of  such  a  shear  at  the  t(^  of  tiie  tower  can  be  handled 
readily  in  the  Rayleigh  solution  by  the  boundary  condition  yf"*  =  at 
z  s  1;  this  boundary  condition  would  be  used  instead  of  ly"'  =  0. 


Base  Shear  and  Overturning  Moment 


There  are  two  static  loads  that  occur  on  a  cantilever  as  a  result  of  earth¬ 
quake  motions.  One  is  the  base  shear,  and  the  other  is  the  overturning 
moment.  These  loads  are 
shown  on  the  distributed- 
mass  structure  of  Rguic  28. 

These  two  loads  are  com¬ 
puted  by  simple  statics, 
and  the  proc^ure  is  pre¬ 
sented  in  this  section. 


In  the  Rayleigh  solu¬ 
tion  for  a  distributed 
mass,  the  acceleration  of 
the  mass  m*  at  the  refer¬ 
ence  level  r  is  found  by 
multiplying  the  spectral 
acceleration  times  the 
normalization  ratio 
Lo/m*.  For  any  given 
shape  function,  the  accel¬ 
eration  of  the  structural 
mass  at  all  other  points 
along  the  height  of  the 
tower  is  found  by  multi¬ 
plying  the  acceleration  at 
the  reference  level  r  by  Figure  28.  Base  shear  and  overturning 
the  shiqie  function.  The  moment 
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molt  tfiea  beconet  die  aeceUradon  fimetion  (Table  3).  Such  a  calcula- 
tioa  for  die  itep-tapered  tower  of  the  preceding  example  is  demonstrated 
in  the  fidlowing  procedure. 


TabtoS 

Aeoalaratlon  Function  Coaffictants 

NonaUbaUM  IMIe  1 J70 

tlniinalliaMnn  natin  l  im 

NenuaanUoa  RaUo  0J0S4 

thoM  t..  ihMa  nrat 
Miap*AnoUoii(TaM»4) 

anain  Swollen  (Tafela  4) 

ttaiaa  Bm  Uawa  ThM 

Btiapa  meUon  (TaMa  4 

a -2.15045 

a -4.03575 

a-7.a0370 

1  Fkatal9p,0<xsi/S  | 

A,-*0.45222 

A, --40.24400 

A, -40.12243 

B,- *0.60004 

B,- -0.25071 

B,- -0.12373 

C,- -0.45222 

C, --0.24400 

C,  - -0.12243 

0,  -  ■fO.00084 

0,  -  40.25671 

0,-40.12373 

1  Saoond  sMp,  1/3  s  X 1 2/3 

A,--»0.40200 

A, -40.27505 

A, -40.14355 

B,- -0.07071 

B,--0.30775 

B,- -0.14300 

C,- -0.75100 

C,- -0.20570 

C,  -  40.12510 

0,-^.05724 

D,-40J20104 

0,  -  -0.12530 

1  Third  atap,  2/3  s  x  S 1 

A,  -  40.44122 

A, -40.32240 

A, -40.17500 

B,- -1.04107 

B,- -0.30010 

B,  -  -0.17326 

C,- -1.24003 

0,-41.34404 

C,- -2.98413 

D,- -1-1.42200 

0,- -1.35001 

0,-4^  3428 

With  these  coefficients,  the  accelerations  are  given  by  the  following 
sets  of  three  discontinuous  equations,  now  termed  the  acceleration 
functions. 

First  acceleration  function 


OSti  t/3:  v  «  0.45222  tin  (2.150450  -  0.66994cot  (2.l5045z) 

-  0.45222 finh  (2.l5045t)  -t-  0.66994 coth  (2.150450 


54 


Chiplnr3  OMributwHUaM Syatsms 


-  ATSiiNM  aimw 

3/3  ittSU  l.MI<T?cflaaiJMW 

.  IJ«N»iWiatS04SO-f  1.4Xlitoaikai304St> 

Secwid  accelenttioB  ftiactkw 

9*t*tAt  V  •  0.a44Mili  (4J337j«  «  0.2S67I cm (4.taS7S0 

>  0L3449tMi  (4.nS75^  025471  each  (4.t357S{) 

l/3ScS2/3e  v- 0^27545  ita  (4.135754 -0Ja77Secc(4J3S75d 

>  029574  alah  (4.435754  029104  eoth  (4.435754 

2/3  S  <  4  I:  V  -  032a494i  (4J35754  -  0349t4GM  (4.43S7S4 

*  I J449444I  (4.435754  >  1.35041  cedi  (4.435754 

Third  acceleratioa  function 

0  4  <4  1/3:  V  -  012243 dc  (7J53744  -  OI2373eM  (7.453744 

-  012243  ddi  (7.453744  *  012373  ccdi  (7.853744 

t/3  4  t  4  2/3:  V  •  OI4355dc  (7.453744  -  0.14399cm  (7.453744 

*  012514  dak  a453744  -  0.12530  codi  (7.453744 

2/3  4  <  4  1:  V  *  0.17344da  (7J53744  -  O.I73:Mcm  (7.453744 

>  294414  dah  (^.$53^160  *  2.94424  CMb  (7.453744 

Rfore  29  io  o  graph  of  the  acceleratitm  functions  along  with  their 
SRSS  condiined  accelerations.  The  CQC  method  of  combining  sh^  fiiiw* 
dons  d^nerates  to  the  SRSS  method  in  all  of  the  shape  functions  encoun* 
tered  so  far.  The  ratio  Tj^*]  is  so  small  that  there  is  never  any 
cross>interaction. 


Figure  29.  Aeoeleration  functions 


CtwiNws 
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The  acceleratioo  is  now  coopted  each  10  ft  along  the  height  of  the 
tower  for  each  of  the  three  acceleration  fiinctiont.  The  calculations  are 
made  from  the  foregoing  equations,  again  using  the  mathematics  software. 
The  results  of  these  calculations  are  shown  in  the  following  tabulation. 


The  tabulated  calculations  give  the  three  accelerations  that  occur  at  the 
mithile  of  each  10-ft  increment.  Each  of  these  accelerations  is  the  maxi¬ 
mum  spectral  acceleration  for  its  particular  shape  function.  The  three 
maximum  accelerations  at  each  increment  are  then  combined  into  a  resul¬ 
tant  acceleration  using  SRSS.  The  force  at  the  middle  of  each  10-ft  incre¬ 
ment  is  the  mass  of  the  10-ft  segment  times  the  combined  acceleration  at 
that  level. 


56 


ChaptwS  OMrituM-Mass  System 


TMourtuevn 

fpwiwitjt  mf  hv9t  im  tefwwrcii  be  readily  calcaltiDd  from  these 
faww—d  h  hKThwM  at  tfw  tort  c^anai  ia  A>  triwriitkML  Figure  30 
dwwt  •doeieli  of  the  renlie. 


KmcEMMra 


SHEAR  OMOflAM 


MOMairOMORAM 


Figure  30.  Forces  and  moments  on  the  step-tapered  tower 


The  renuiiiing  calculation  to  be  made  is  that  for  the  displacement  of 
the  tower  after  it  has  sustained  some  cracking  firom  the  maximum  design 
earthquake  (MDE).  The  generalized  mass  m*  and  the  normalization  fac¬ 
tor  Lg  are  unchanged  in  this  event,  but  since  El  is  assumed  to  be  reduced 
by  half,  die  stiffness  k*  is  reduced  by  l/>^.  ^th  these  revisions,  the 
Rayleigh  constants  and  the  natural  frequencies  become  those  listed  in  the 
following  tabulatirm. 


h  Constuits  for  MDE  Dioptacemonts 


Flt9t  ShBp#  funotton 

SeceiNl  Sbap*  Rmetien 

Third  Shaps  Funetien 

HI*  >  90,560  skigi 

m*.  77,600  akigi 

m*  -  79,700  slugs 

k*»25.6x10*Kim 

k*-602x10*lim 

k*  -  4,325  X 10*  bm 

L.«  151,215  skiOi 

L„a -79,460  slugs 

L„  -  50,540  slugs 

LnAn*- 1.670 

IVHl*- -1.022 

LVm*  «  0.634 

««1641  radteM 

«>66.0rad/ssc 

«-233rad/Sse 

T  •  0 J74  see 

T- 0.071  ssc 

T  -  0.027  ssc 

ClH|Nw3  tUsaaiiilid  IKesi  Oysisnw 


The  revised  spectral  accelerations  are  found  from  Figure  17  for  the 
periods  just  calculated. 


past  Shape  PuiMlioiit  q 

Saeend  Shape  Function,  a 

Third  Shape  Function,  g 

8^-1.128 

8a -0.770 

8a -0.572 

^th  these  revised  values  of  spectral  acceleration,  the  entire  procedure 
introduced  earlier  is  repeated  to  find  the  revised  equations  for  the  three  ac¬ 
celeration  functions,  llie  new  accelerations  are  then  recomputed  ^m  the 
revised  acceleration  functions  at  10-ft  intervals  along  the  tower.  The  dis¬ 
placements  corresponding  to  these  new  accelerations  are  computed  simply 
by  dividing  the  accelerations  by  oo^.  The  results  of  such  calculations  are 
given  in  the  following  tabulation. 


III 

Maplaoemants  Under  MDE,  ft 

Diaplacanienta 
Under  QBE,  ft 

Fhat 

Second 

Third 

suss 

Uncrackod 

ft 

Function 

Function 

Function 

ComMnod 

Section 

178 

•fO.2058 

-0.0029 

♦0.0002 

0.2055 

0.1455 

165 

f 0.1 884 

•0.0020 

♦0.0001 

0.1884 

0.1334 

155 

♦0.1713 

•0.0012 

•0.0000 

0.1713 

0.1213 

145 

♦0.1544 

•0.0005 

•0.0001 

0.1544 

0.1097 

135 

♦0.1377 

♦0.0002 

•0.0001 

0.1377 

0.0975 

125 

♦0.1214 

♦0.0008 

-0.0001 

0.1214 

0.0859 

115 

♦0.1056 

♦0.0013 

-0.0001 

0.1056 

0.0747 

105 

♦0.0904 

♦0.0016 

-0.0001 

0.0904 

0.0640 

95 

♦0.0758 

♦0.0018 

•0.0000 

0.0758 

0.0537 

85 

♦0.0821 

♦0.0019 

♦0.0000 

0.0821 

0.0440 

75 

♦0.0495 

♦0.0018 

♦0.0001 

0.0495 

0.0350 

65 

♦0.0381 

♦0.0018 

+0.0001 

0.0381 

0.0270 

55 

♦0.0281 

♦0.0014 

♦0.0001 

0.0281 

0.0200 

45 

♦0.0194 

♦0.0010 

♦O.0QO1 

0.0194 

0.0137 

35 

♦0.0121 

♦0.0007 

♦0.000 1 

0.0121 

0.0086 

25 

♦0.0064 

♦0.0004 

♦0.0001 

0.0064 

0.0044 

15 

♦0.0024 

♦0.0002 

♦0.0000 

0.0024 

0.0017 

5 

♦0.0003 

♦0.0000 

♦0.0000 

0.0003 

0.0002 

0 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

58 


Chapters  Oistrlbuted-Mass Systems 


Noli  thtt  is  dw  tabalirion.  dw  dme  revised  eoeelerelioB  feacdoas 
kaw  baea  divkM  by  md  have  duu  beconee  three  diiplaceBieat  fiuic- 
tkau.  The  duee  di^daceamts  at  each  eievadoa  are  con^aed  usiag 
SRSS  aa  showa. 

The  displaceaienta  fw  this  particular  tower  under  the  MDE  are  com¬ 
pletely  dominated  by  the  first  displacement  function.  For  the  sake  of  com¬ 
parison,  the  displacements  under  the  OBE  were  calculated  separately  and 
are  listed  in  die  last  column  of  the  tabulation.  The  displacements  under 
the  OBE  were  similarly  found  to  be  completely  dominated  by  the  first 
displacement  function. 
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4  Multistory  Lumped-Mass 
Systems 


Multiple  Lumped-Mass  Systems 

The  Rayleigh  solution  piesented  in  the  previous  sections  applies  to  sys¬ 
tems  having  a  distributed  mass.  The  mass  in  those  systems  may  be  uni¬ 
formly  distributed  as  shown  in  Figure  31a;  it  may  vary  in  steps  as  shown 
in  figure  31b;  or  it  may  vary  continuously  as  shown  in  Figure  31c.  While 
the  computation  of  the  generalized  mass  and  the  generalized  stiffness  may 
be  quite  complex,  the  use  of  mathematics  software  to  perform  the  integra¬ 
tions  makes  such  solutions  readily  possible  in  modem  practice. 

However,  not  all  systems  are  distributed-mass  systems.  Rigid  frame 
buildings,  for  example,  have  the  building  masses  lumped  at  each  floor, 
with  comparatively  slender  columns  and  negligible  masses  between  the 
floors.  Such  a  structure  is  shown  schematically  in  Figure  32.  When  the 
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Figure  32.  Lumped-mass  system 


building  is  deformed  laterally,  as  in  an  earthquake,  all  of  the  restoring 
force  is  provided  by  flexure  in  the  columns. 

The  multistory  lumped-mass  system  shown  in  Figure  32  is  not  a  single- 
degree-of-freedom  system.  Later  discussions  will  show  that  this  three- 
mass  system  will  have  three  degrees  of  freedom,  three  natural  frequencies, 
and  thiM  defmmation  shiqies.  It  is  classed  as  a  multiple-degree-of-lteedom 
(MDOP)  system. 

As  previously  stated,  however,  the  distributed-mass  systems  treated  in 
the  Rayleigh  solution  were  also  MDOF  systems,  but  it  was  possible  to 
examine  their  dynamic  responses  by  considering  only  one  degree  of  free¬ 
dom  (one  shape  function)  at  a  time.  Therefore,  the  multiple  lumped-mass 
systems  may  also  be  examined  by  considering  only  one  degree  of  freedom 
(one  mode  shape)  at  a  time.  The  shape  of  the  deformation  curve  in  a 
lumped-mass  system  is  called  the  mode  shape  and  is  analogous  to  the 
shape  function  of  the  Rayleigh  solution. 

In  most  existing  building  systems,  the  floor  beams  are  much  stiffer  in 
flexure  than  the  columns.  The  structure  will  therefore  deform  essentially 
as  if  the  columns  were  fixed  at  each  floor  and  can  be  analyzed  as  the 
close-coupled  system  shown  in  Figure  32  with  a  deformed  shape  as  shown 
in  Figure  33.  Tlie  structural  idealization  is  designated  as  close-coupled  be¬ 
cause  only  the  lateral  displacements  at  each  floor  level  are  allowed  to  in¬ 
fluence  tte  distribution  of  inertial  forces  at  the  various  floor  levels.  The 
current  seismic  design  practice,  however,  is  to  provide  columns  that  are 
stronger  than  the  beams.  The  stronger  column  system  prevents  early  de¬ 
velopment  of  plastic  hinging  in  the  columns.  Plastic  hinging  in  columns 
can  lead  to  collapse  of  the  building.  When  analyzing  strong  column/weak 
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Figure  33.  Lateral  deformations  of  a  lumped-mass  system 


beam  systems  for  earthquake  ground  motions,  it  is  desirable  to  use  a  far- 
coupled  analysis  procedure.  A  far-coupled  structural  idealization  is  one 
that  includes  the  effects  of  beam/colunm  joint  rotations  on  the  distribution 
of  inerial  forces  at  the  various  floor  levels.  The  close-coupled  system, 
however,  is  used  in  the  following  example  to  analyze  lumped-mass  build¬ 
ing  systems  for  earthquake  ground  motions. 

There  are  two  cases  of  deformation  that  occur  in  the  columns  shown  in 
Figure  33,  that  of  a  beam  fixed  at  both  ends  and  that  of  a  beam  fixed  at 
only  one  end.  The  two  cases  are  shown  in  Figure  34  along  with  the  formu¬ 
las  for  calculating  the  deflections  produced  by  a  transverse  load  P  applied 
at  the  ends.  The  beams  depicted  in  Figure  34  may  look  more  familiar  if 
they  are  turned  90  deg  counterclockwise. 


For  the  beam  loaded  at  one  end  only,  the  spring  constant  k  is  computed 
as  the  force  per  unit  deflection. 


k 


P 

A 


3EI 


Ib/in. 


(60) 
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Figure  34.  Force-deflection  relations  in  beams 


Similariy  for  the  beam  fixed  at  both  ends,  the  spring  crastant  k  is  com- 
pntedas 

.  p  mi (61 

ifc  =  T  =*  — r*  lo/in. 

A  i} 

The  spring  constant  k  is  conunonly  called  the  stiffness  of  a  member  in 
flexure,  whether  a  beam  or  a  column.  It  is  the  restoring  force  per  unit  de¬ 
flection  that  will  be  exerted  by  a  beam  or  a  column  having  a  flexural  mo¬ 
ment  of  inertia  I. 


Modal  Analysis  of  a  Lumped-Mass  System 

Figure  35  shows  a  frame  of  a  typical  MDOF  lumped-mass  building. 
The  stiffhess  k  shown  between  two  lumped  masses  is  the  sum  of  the 
stiffnesses  of  all  columns  that  undergo  the  same  lateral  displacement.  At 
each  level  i.  the  mass  is  denoted  m],  the  displacement  of  the  mass  is  de¬ 
noted  y],  and  the  forcing  function  on  the  mass  is  denoted  F{(t). 

Forces  are  summed  on  the  undamped  flee  body  in  Figure  35  as  follows: 

-  *<-!./  -  ^r-i)  +  *u+i  ^^1+1  -  ® 

When  the  forcing  function  does  not  exist,  the  resulting  mode  shape  at 
any  level  i  will  be  that  of  free  viluation  of  an  undamped  system  supported 
laterally  by  the  masses  above  and  below  it.  The  equation  is: 
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Figure  35.  Forces  acting  on  a  lumped  mass 


*  0  (63) 

The  solaticm  to  Equation  53  is.  as  before: 

y^  *  Of  sin  (cor  -  a)  (64) 

where  aj  is  the  maximum  amplitutte  of  the  displacement  For  this  value 
ofy,: 

=  -Oj  CO*  sin  (car  -  a)  =  -co^i  (®5) 


These  values  for  y|  and  are  substituted  into  Equation  S3  to  find,  with 
Fl(t)-0: 

i  (a,  -  a^.,)  +  (a^^,  -  ap  =  0  < 

In  a  more  easily  remembered  form,  the  equation  may  be  written: 


^btlow^btlow  ^above^above  ”*  ® 


(67) 


where 
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aito  ht  cailbd  die  ‘‘ttwe-iaaplitade  equation."  Combined  appropriate 


•oftwaie  for  aolviaf  siaiiittraeoaa  eqaationa,  its  use  permits  a  very  simple 
and  direct  ablation  for  the  aataral  fr^oencies  and  mode  shapes  of  lumped- 


mass  systems.  The  following  exan4>le  illustrates  its  use. 


Example  Solution  for  a  Three-Story  Building 

Given:  Three-story,  lumped-mass  system  as  shown  below. 

Find:  1)  Natural  periods  of  oscillation,  and 

2)  Mode  shapes  corresponding  to  each  period 


Solution: 

Masses  and  stiffiiesses  are  computed  for  each  level. 

mi  s  490,000/32.2x12  »  1,268  Ib/in/sec/sec 
m2  «  435.000/32.2x12  s  1.126  Ib/in/sec/sec 
m3  s  34S,00a/32.2x  12  =  893  Ib/in/sec/sec 

*QJ  -  iO/t?  m  3  X  30 X  10^  (SOS  +  9S0 •••  9S0 S09/(I6 X  12)^  -  ytJOOiU/in 
*12  -  -12x30xl(/^(220-t-400-f4004  22(V(I2  x  12)^  -  I49.49S  Ib/ln 

*23  -  VtEI/1?  -  12  X  30  X  I«f^  (160  +  300  300  +  16(V(I2  x  12)^  >  110.9180/111 
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The  coiqpirted  values  afe  diowu  cm  a  lomped-mass  skiMch. 


K>aMN*IJM 


The  three-force  equatioas  are  written  for  each  mass. 

Mass  No.  1:0  4-  (1.268«s2  .  i85,000)a,  +  t49.000a2  «  0 

Mass  No.  2: 149.000aj  4-  (1.126<d^  -  260.000)a2  4^  1 1  l.OOOsj  *  0 

MassNo.3:  111.000s2  4-(  893cf|2- ni.000)a3  4-  0  -0 

The  result  at  this  point  is  three  equations  in  four  unknowns,  with  the 
variable  appearing  as  a  coefRcient.  As  in  the  distributed-mass  solu¬ 
tion,  the  lunp^mass  solution  is  an  eigenfunction,  having  one  eigenvalue 
oftrfor  each  equation.  The  amplitude  a3  at  the  top  level  of  the  building 
is  now  set  equal  to  unity,  reducing  the  set  of  equations  to  three  equations 
in  three  unknowns.  All  of  the  other  anq>Iitodes  will  then  be  found  as  mul- 
rifdes  of  amplitude  a3,  termed  normalized  amplitudes. 

The  three  simultaneous  equations  are  solved  using  conventional  mathe¬ 
matics  software.  No  matter  how  the  set  of  equations  is  solved,  the  solu¬ 
tion  for  <n  will  always  be  an  n^  degree  polynomial  in  Thus,  there  will 
always  be  n  solutions  fmr  the  set  of  n  equations.  This  set  of  equations  will 
therefore  have  three  solutions,  one  for  each  natural  frequency  m.  The 
starting  value  of  to  (required  by  the  software)  is  varied  manually  until  all 
three  solutions  are  fou^.  (See  Appendix  A  for  an  alternate  solution  proce¬ 
dure  for  <D  using  tnatrix  methods.) 


•i «  3.1SB  radteae 

■  11 J02  ladteao 

-  IBJOS  radteac 

T  -  IJSesae 

T  «  0.842aao 

T  -  0.332aac 

a,  -  40.791 

at-<0J0B 

a, -1-1.046 

•fm40MO 

a,- -0.083 

a2--l.87S 

a,o-»1.000 

a,  « -f-1.000 

a,-i-l.00C 
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A  lignificait  featuie  leveaied  by  Ihu  eolation  is  that  the  lateral  deflec- 
tioo  of  tfM  fint-floor  mass  is  quite  high.  Almost  80  percent  of  the  total  de- 
flectkm  in  die  first  mode  shape  ocean  in  the  fint  floor.  One  cause  of  this 
high  defleetkm  is  that  the  first-floor  columns  are  longer  than  the  other  col- 
nmns,  a  oommoB  design  fhataie.  Fhither,  die  fint-floorcohunns  are  hinged 
St  their  bases  rather  than  being  fixed,  making  them  much  more  flexibte. 

As  a  point  of  interest,  if  the  colunms  at  the  fint  level  could  be  fixed  at 
their  bases,  the  boilding  would  be  much  more  rigid,  and  the  fint-floor  de¬ 
flections  would  be  reduced.  The  following  calculations  show  why  this  is 
so. 

fixed-etdumn  bases  at  ground  level,  the  stiffness  factor  at  the  fint 
level  becomes: 

kQx  •  \2EI/L^  «  148,000 /^/u. 

This  value  is  used  in  the  three-force  equation  for  mass  No.l, 

Mass  No.  1:0+  (1.268<d2  -  297.000)a,  +  149.000a2  «  0 

With  this  change  in  the  set  of  three-force  equations,  the  computer  solution 
now  yields: 


The  foUowing  dcetches  show  the  node  thapet  for  diit  more  rigid  stracture. 


«M(  MM  MM 


\inien  this  second  (OKm  rigid)  response  is  compared  with  the  first 
(moire  flexible)  re^KMise,  it  can  be  seen  that  the  deflections  of  the  flrst 
mass  are  reduced  considerably  in  the  first  two  modes,  with  but  little 
change  in  the  third  mode  deflection.  In  the  flrst  mode,  the  deflection  is  re¬ 
duced  to  about  half  of  what  it  was.  Assuming  that  all  the  upper  floors  de¬ 
flect  the  same  amount  in  both  solutions,  the  total  lateral  deflection  of  the 
building  has  been  reduced  approximately  35  percent.  The  most  striking 
difference,  however,  is  in  the  two  periods.  The  first-mode  period  of  the 
more  flexible  building  is  almost  67  percent  longer  than  the  period  of  the 
more  rigid  building.  This  period  elongation  brings  the  response  of  the 
more  flexible  structure  well  below  the  peak  values  of  spectral  accelera- 
ticms  given  by  the  ATC  3-06  corves  shown  in  Figure  17.  Since  a  further 
elongation  of  the  period  would  reduce  the  spectral  accelerations  even 
more,  such  an  alternative  seems  worth  pursuing. 

The  results  of  the  fwegoing  comparisons  suggest  that  the  design  might 
be  in^MPOved  by  base  isolation.  For  that  case,  the  ground  floor  would  be 
raised  off  the  foundations  at  each  column  by  isolators,  each  having  a  stiff¬ 
ness  of  5,000  Ibfln.  Such  a  base  isolation  system  is  shown  in  the  follow¬ 
ing  sketch.  For  best  utilization  of  the  base  isolation,  the  ground-floor 
columns  me  fixed  (rather  than  hinged)  to  the  ground-floor  girders. 
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(lloarltfhN»atS40g0Q0lMi.  TlM»^ovB,dM 
rttl^QOftMM/Mc.  ThtstifAMttorewhoftlMfMr 
iMlitan  ift  aet  •!  3,000  AAo,  for  t  total  itiffiMaf  at  tlie  isolittor  level  of 
20,000  Ih^.  The  tltt«e>foice  eqiiatioas  become: 


Mast  No.  0:  0  •••  (l,400«s2  •  168.000)ao  148.000a,  -  0 

Mass  No.  1: 148.0000^  -i-  (1.267ci>2  -  297.000)a,  +  149.000a2  »  0 

Mass  No.  2: 149.000a,  (I.125m2  >  260.000)82  ^  1 1  l.OOOa^  «  0 

MassNo.3:  lll.OOOajX  892(0^  -  111.000)83  •(-  0  -0 

The  solutions  of  these  four  equations  yield  the  following  results  for  the 
diree>8tory  lumped-mass  building  with  base  isolators. 


•,«  i.9aeradnM 

•^>  S.S37rad/SM 

>  15.036  rad/aae 

T  •  3.104  MO 

T  -  0.711  MC 

T  -  0.416  aae 

a,. -0.797 

ao- 40.771  . 

a,  -  40.917 

a,  >  -0.316 

a,  >  -0.774 

1^.40 J69 

1^.40.372 

a^  >  -0.819 

a,  >41.000 

a|  >41.000 

a, -41.000 

The  following  sketches  show  the  shape  functions  with  base  isolators, 
with  ground-floor  columns  flxed  at  their  bases. 
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The  base  isolators  take  about  80  percent  of  the  total  lateral  displace¬ 
ment,  thefd>y  sharply  decreasing  the  relative  displacements  in  the  struc¬ 
ture  shown  in  the  previous  sketch.  Further,  the  isolators  signiticantly 
increase  the  period  T,  bringing  the  spectral  resp<mse  from  the  ATC  3-06 
curves  well  down  from  the  peak  value. 
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Thus,  it  can  be  concluded  that  a  principal  advantage  of  base  isolators  is 
that  they  lengthen  the  natural  period,  thereby  reducing  the  spectral  acceler¬ 
ations  and  consequent  inertial  forces.  One  of  the  disadvantages  is  that 
large  overall  displacements  must  be  accommodated  in  the  building  design. 


Modal  Participation  Factor  for  a  Lumped-Mass 
System 

The  n-story  building  riiown  in  Figure  36  is  subject  to  earthquake  excita¬ 
tion.  The  equations  of  motion  for  the  building  are  obtained  by  summing 
forces  on  the  free  body  of  each  lumped  mass. 


Figure  36.  Typical  lumped-mass  building 


The  equations  of  motion  are: 

"*1^1  *01  0^1  -  *12  0^2  -^i)  =  ® 

+  *12  ^2  ■  *23  ~^2)  =  ® 

.  (68) 

+  *«-2.«-l  (yn-l  -yn-2^  "  *n-l.«  O',  =  0 

As  was  done  with  SDOP  systems.  Equation  6,  the  coordinate  displace¬ 
ment  y]  is  transformed  into  the  displacement  relative  to  the  ground  by  the 
transformation: 


70 


CtMpIsrS  MuWstOfy  Lumpxt-Msw  Systsim 


f  ' 


l! 


Sriwtotkwi  into  dM  e<pwtkMM  of  modoa  yields: 


•A  ♦  \ii»i  -  »u 

■V2  +  *12  -  *23  -  -2*J» 

-  (70) 


•%_  iS,>  1  ♦  1  ‘"e-  I-*ii-2>  -  *«- 1.*  <%-•«- 1>  •  ^11- 1“«« 

%%♦*»-  l.ii  ‘•n  '  ^1- 1>  - 


The  system  of  eqnetioiu  is  nncottpled  by  transforming  the  diH>lace- 
ments  oi  into  a  set  of  time  ftinctioos  Z|<t)  times  the  undamped  nu^s,  one 
ftmctioa  Z|  for  each  mode: 


“l  -  «il«i  +  *124  ♦  «I3«3  + . 

-2  -  «2I*1  +  «2A  +  «23*3  +  . 

.  (71) 


“«  *  "el^I  +  «a2*2  +  "1.3^  + 


where  the  second  subscript  for  the  amplitude  a  denotes  the  mode. 


These  emotions  are  substituted  back  into  the  equations  of  motion. 

^di  some  rather  complex  matrix  numipulations  (^  1991),  the  final  solu¬ 
tion  for  the  time  functions  Z{  for  each  niode  are: 

z,  +  oft,  =  r,u^(r) 

.  (72) 

'if  +  of  Zt»  Tftjjt) 

The  term  Fj  for  each  mode  i  is  given  by: 


H  H 

Where  -  X  •«*  «  Z  w/pj 

/-I  /-I 


(73) 
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The  symbol  Oj,  denotes  the  unforced  amplitude  of  the  mass  oscillating 
inthe>modil 

Hie  right-hand  sides  of  Equation  72  are  the  same  at  all  levels  of  the 
building  for  any  given  mode.  The  factor  may  therefore  be  viewed  as  a 
multiplier  on  the  ground  motions  dgO).  Consequently,  when  ground  mo¬ 
tions  are  substituted  into  the  transformation  equations,  Fj  also  becomes 
the  multiplier  on  the  relative  amplitudes  at,  a2....  a^.  The  factor  Fj  is  the 
multiplier  that  tranforms  the  relative  values  of  modal  displacements  into 
absolute  values. 

niysically,  the  multiplier  Fj  indicates  the  increase  or  decrease  of  the 
ground  motion  that  occurs  when  the  various  modes  of  the  time  functions 
zj  are  summed,  each  of  which  may  be  positive  or  negative  at  any  given  in¬ 
stant.  F}  is  called  the  modal  participation  factor,  and  is  directly  analo¬ 
gous  to  Ae  integral  forms  developed  in  the  Rayleigh  method,  right-hand 
side  of  Equation  52. 

The  absolute  values  of  displacement,  velocity,  and  acceleration  for 
each  mass  can  now  be  expressed  in  terms  of  the  spectral  values  of  Sq,  Sy, 
and  for  each  mode  shape  i. 

u..  =  F.ffl.-So- 

jt  mu  I  ^jt  Dj 

s  F.®..Sw.  (74) 

*7»  mu  » “V;  ^  ' 

ii..  aF.ffl.S. 

ji  m$x  t^jt  A] 

The  values  of  Uj,  Uj,  and  (ij  for  the  various  modes  can  be  combined 
either  by  SRSS  or  by  CQC  to  obtain  the  Hnal  combined  values.  If  SRSS 
is  used,  the  combined  values  are  given  by: 

«  =  Vs 

jmax  ^ 

it.  =  Vs 

jmax  M 

y=i 

The  force  fj  acting  at  any  level  j  is  simply  the  acceleration  at  that  level 
times  the  mass  at  that  level: 

fj  =  ^'hmax  (76) 
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Chapter  4  Multistory  Lumped-Mam  Systems 


Tlie  overturaiag  momeot  on  a  lomped-maas  lyttem  can  be  found  simi¬ 
larly.  The  force  ft  acting  at  any  level  has  already  been  found  and  defined 
by  Equation  39.  ilw  overturning  moment  Mg  about  the  base  of  the  struc¬ 
ture  caused  by  the  forces  fj  is  the  sum  all  these  forces  fj  times  their 
height  yj  idiove  the  base: 

"  (78) 

f-1 


Equations  77  and  78  yield  the  design  values  for  base  shear  and  over¬ 
turning  mommit  fw  a  Innqwd-mass  system  subjected  to  earthquake  load¬ 
ing.  Prom  this  point  onward,  the  remainder  of  the  analysis  is  one  of 
tmiinaty  statics.  The  force  at  any  level  y  in  the  structure  is  given  by  Equa¬ 
tion  76.  The  shear  at  any  level  y  is  then  simply  the  static  sum  of  forces 
above  that  level.  The  moment  at  any  level  y  is  the  static  sum  of  moments 
due  to  the  forces  above  that  level. 
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Appendix  A 

Matrix  Soiution  for  a  Three-Mass 
MDOF  System 


The  general  equation  of  motion  fw  a  lttmped*mass  system  is  given  by 
Equation  66  in  the  main  text  When  there  is  no  forcing  function,  the  equa¬ 
tion  becomes: 


+  ('"X  -  X*/)  "i+i  =  ® 

In  alternate  f<ma. 


When  applied  to  a  lumped-mass  structure  such  as  that  shown  in  Hguie  A-1, 
the  set  of  equations  becomes: 


0  —  (kQi  +  ^12)^1  *12®2  ” 

*13^1  ~  (*12  *23^  ” 

*23^2  ”  ^*23  *34)  ®3  *3  A  ~  ”“*^''*3^ 

.  (A-3) 


In  matrix  notaticm,  the  set  of  equations  is  given  as: 

[jr  -  a^M]m  =  0 

where  K  is  the  stiffness  matrix.  M  is  tiie  mass  matrix,  and  a  is  the  ampli¬ 
tude  matrix. 


The  nontrivial  solution  for  Equation  A-4  (i.e.,  fw  a  #  0)  requires  that 
the  determinant  of  die  coefficient  matrix  most  be  zero: 


AoiimaicA  lirei>  SoMlBn  for  a  Thna4laM  MDOF  SnlMi 
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Figure  A-1.  General  iumped-maee  building 


K  -  *  0  (A-5) 

Equation  A>5  defines  an  eigenproblem.  The  solution  of  the  determinant 
yields  a  polynomial  equation  of  the  n^  degree  in  a>^.  The  n  roots  of  tel^ 
are  called  the  eigenvalues  of  the  problem.  A  different  set  of  values  for  the 
ampUtudes  a^,  a2***'^  occur  for  each  eigenvalue. 

As  an  example,  the  solution  will  be  applied  to  the  three-mass  structure 
solved  earlier  in  Chipter  4  of  the  main  text  and  shown  again  in  Figure  A-2. 
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AcpwvaxA  MaMx  SoMlon  fw  a  TlviM44aM  MDOF  SytlMi 


Tie  ttUbeu  netrix  K  Ik  given  by 


-186,000  •(•149.000  -i-  0 

Km  149.000  -260.000  -1-111.000 

0  -(-111.000  -111.000 

The  mask  nutrix  M  ik  giv«i  by 

-1.268  0  0 

Mm  0  -1.126  0 

0  0  -893 

The  coefficient  matrix  ik  found  from  these  matrices: 

■lS6.00l>fiaSS«^  ••■149^  0 

i49jaM  ^tfaooofi.i26«^  -t-nixoo 

0  +lll/»0  -111.0004- S93»^ 

The  determinant  is  expanded  by  minors  to  find: 

(m*)^  -  502(«^)*  +  53,090(oo?)  -  481,130  =  0 
The  roots  of  the  cubic  equation  are  found  by  trial  and  error: 

<0|  =  ±3.16  rad/sec  ;  0012  =  ±11.61  rad/sec  ;  (Uj  =  ±18.91  rad/sec 

For  the  first  mode,  <D  =  3.16  rad/sec  and  the  matrix  becomes: 


r-173.350  -1-149.000  0  ] 

[J:  -  ooi^Mla  =  149,000  -248.770  -i-lll.OOO  (*2 

0  -mi.OOO  -102,090j  oj 

The  set  of  equations  is  solved  simultaneously  to  find: 

Uj  =  -(0.791  Uj  °  -^-920  =  -(-1.000 

The  scdution  is  repeated  for  the  second  and  third  eigenvalues  to  find 
the  sectmd  and  diitd  mode  shiqies.  The  solutions  are  summarized  below: 

mj  =  3.16  rad/sec  otj**  11-61  rad/sec  (d^  =  18.91  rad/sec 

a|  =-(0.791  a}«-0.80S  a.  =-(1.045 

•2  -•(0.920  a2  = -0.083 

aj  =-(1.000  a3  = -(1.000  ajs -(1.000 

These  solutkms  are,  of  course,  identical  with  those  obtained  for  this 
pfoblmn  in  Oe  exanqrte  given  in  Charter  4  of  the  main  text 
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Appendix  B 
Notation 


A,v,D 
A,  B,C,  AD 
A* 
c 
d 
O' 

E 

El(x) 

f 

fj^dt 

fi^mi 

U^ku 

I 


PMk  groond  acceleratira 

Coostantt  in  Chapter  2  (see  Ihble  1  for  values) 

Constants  of  integratitm  in  Chapter  3  (Equation  56) 

Maximum  pseudoacceleration  m^D'  of  the  SDOF  mass 

Damping  constant 

Peak  grmmd  displacement 

Maxinmm  relative  displacement  between  ground  and 
mass  of  an  SDOF  system  having  a  natural  frequency /at 
whatever  point  in  time  it  occurs 

Modulus  of  elasticity 

Flexural  stiffness 

Frequency  of  oscillatimi 

Danq>ingfmce 

Inertial  force 

Restoring  force 

Moment  of  inertia 


k  Stiflhess;  also,  spring  constant 
K  Stifhiess  matrix  . 

KE  Kiamkenrngy 


Lt/m*  DimensioBleM  multiplier  applied  to  ground  acceleration 


m  Mass 

Mass  per  unit  lengdi 
M  Moment  on  die  cross  section 
MDE  Maximum  design  earthquake 
OBE  Operational  basis  earthquake 
jy(t)  Externally  tqiplied  dynamic  force 
p(t)  Harmonic  force 

Cross'funcdon  coefficient  in  which  i  and  j  are  any  two 
shape  functions 

P0  Maximum  magnitude  of  the  forcing  function 
P(t)  lime-dependent  force  of  magnitude  -  m  • 

PE  Stored  potential  energy 
PGA  Peak  ground  acceleration 
PGD  Peak  ground  displacement 
PGV  Peak  ground  velocity 
r  Reference  level 

SA  Absolute  accelermion  response  spectrum 
5^  Spectral  pseudoacceleration 
Sjj  Spectral  maximum  displacement 
Sy  Spectral  maximum  pseudovelocity 
t  Hme 

T  Undamped  period  of  vibration 
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Tq  Damped  period  of  vibration 
T  undamped 

u(t)  Harmonic  displacement 


Appsndx  B  NotsSon 


m^t)  BaitfKpuke  froniid  motion 
Maximum  di^lacement 
Ug  Initial  displacement 

Ugf  Diqdaonnent  that  would  produce  if  ai^died  as  a  static 
force 

ut(t)  Total  diq>lacement  of  mass  relative  to  its  at-rest  position 
ito  Initial  velocity 
ujit)  Total  acceleration 

V  Peak  ground  velocity 

V  Shear  on  die  cross  section 

V*  Maximum  pseudovelocity  tal/  of  the  SDOF  mass 

P  Constant  equal  to  c/2mm 

P  PracdtHi  of  critical  damping 

r,-  Modal  participatimi  factor  for  mode  i 

0  Rotadcm  of  the  cross  section  about  the  neutral  axis 

9  Lagging  phase  angle  between  motion  of  ground  and 
motion  of  mass 

v(t)  Hme  function 

Vf,  Maximum  displacement  at  reference  level  r 
qr  Displacement  fnnctimi 
ip  Shape  function 
o>  Natural  circular  frequency 
os  Undamped  angular  frequency  of  vibration 
os^  Damped  circular  frequency 
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«£>  aiignlar  freqaeacy  9i  vibration 

ft  Forcing  frequency 
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straint  Processing,  Volumes  I  and  II 

A  Case  Committee  Study  of  Rnite  Element  Artalysis  of  Concrete 
Flat  Slabs 

User’s  Guide:  Computer  Program  for  Two-Dimenshnal  Analysis 
of  U-Frame  Structures  (CUFRAM) 

User’s  Guide:  For  Conaete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-83 

Rnite-Element  Method  Package  for  Solving  Steady-State  Seepage 
Problems 

User’s  Guide:  A  Three  Dimension^d  Stability  Analysis/Design 
Program  (3DSAD)  Module 
Report  1 :  Revision  1 :  General  Geometry 
Report  2:  General  Loads  Module 
Report  6:  Free-Body  Module 


Dide 

Jun1982 

Jan  1983 

Jun1983 

Jul1983 

Sep  1983 
Sep  1983 

Oct  1983 
Jan  1984 

Aug  1984 

Sep  1984 

Sep  1984 

Oct  1984 

Jun  1986 

Jan  1987 
Apr  1987 
May  1987 
May  1987 

Jun  1987 

Jun  1987 
Sep  1989 
Sep  1989 


(Continued) 


WiraiWAYS  iXPEMICIIT  STATION  REIK>^ 
PUBU8HED  UNDER  THE  COMPUTER-AIDED 
STRUCTURAL  ENQINESraNQ  (CASE)  PROJECT 


Instruction  Report  ITL-87-4 
Tschrycal  Ftoport  ITL-87-4 


Instruction  Report  QL-87-1 

Instruction  Report  ITL-87-S 
Instruction  Report  rTL-87-6 

Technical  Report  ITL-87-8 

Instruction  Report  ITL-88-1 

Technical  Report  ITL-88-1 

Technical  R^rt  ITL-88-2 

Instruction  Report  ITL-88-2 

Instrtjction  Report  ITL-88-4 

Instruction  Report  QL-87-1 

Technicai  Report  m.-89-3 
Technical  Report  iTL-89-4 


(Cor^nued) 

Title 


User’s  Guide:  2-0  Frame  Analysis  Link  Program  (LINK2D) 


Finite  Element  Studies  of  a  Horizontally  Framed  Miter  Gate 

Report  1 :  Initial  and  Refined  Rnite  Element  Models  (Phases 
A.  B,  and  C).  Volumes  I  and  II 
Report  2:  Simplified  Frame  Model  (Phase  D) 

Report  3:  Altsmete  Configuratior  Miter  Gate  Rnite  Element 
Studtes-Open  Section 

Report  4'  Alternate  Configuration  Miter  Gate  Rnits  Element 
^  Studies-Closed  Sections 

Report  5:  AHemate  Configuration  Miter  Gate  Rnits  Element 
Studiee-Additionai  Closed  Sections 
Report  8:  Bastic  Budding  of  Girders  in  Horizontally  Framed 
Miter  Gates 

Report  7:  Apfriication  and  Summaiy 


User's  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  I, 
User's  Manual 


Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Criteria  Specifications  for  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 
Gates  (CMITER) 

Procedure  for  Static  Analysis  of  Gravity  Dams  Using  the  Rnite 
Element  Method  -  Phase  la 

User's  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

User’s  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

User's  Guide  for  Design  and  Investigation  of  Horizontally  Framed 
Miter  Gates  (CMITER) 

User's  Guide  for  Revised  Computer  Program  to  Calculafo  Shear, 
Moment  and  Thrust  (CSMT) 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  II, 
Theory 

User’s  Guide:  Pile  Group  Analysis  (CPGA)  Computer  Group 

CBASIN-Structural  Design  of  Saint  Anthony  Falls  Stilling  Basins 
According  to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures;  Computer  Program  X0098 


Date 
Jun  1987 
Aug  1987 


Aug  1987 

Oct  1987 
Dec  1987 

Jan  1988 

Feb  1988 

Apr  1988 

Apr  1988 

Jun  1988 

Sep  1988 

Feb  1989 

Jul  1989 
Aug  1989 


(Continued) 


WATERWAYS  EXPERfMENT  STATION  REPORTS 
PUBUSHED  UNDER  THE  COMPUTER-AIDED 
STRUCTURAL  ENGINEERING  (CASE)  PROJECT 


Technical  Report  ITL-89-5 

Technical  Report  rTL-89-6 
Contrad  Report  ITL-69-1 
Instruclion  Report  rTL-90-1 
Tedmical  Report  m.-90-3 

instruction  Report  rTL-90-6 
instnxrtion  Report  rn.-90-2 
Technical  Report  m-91-3 

Instruction  Report  ITL-91-1 

Instruction  Report  fTL-87-2 
(Revised) 

Technical  Report  rTL-92-2 
Technical  Report  ITL-92-4 
Instruction  Report  ITL-92-3 

Instruction  Report  rn.-92-4 
Instruclion  Report  ITL-92-5 


(Continued) 

Titie 

CCHAN-Structucal  Design  of  Rectangular  Channels  According 
to  Corps  of  Engineers  Criteria  for  Hydraulic 
Stnjctures;  Computer  Program  X0097 

The  Response-Spectrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Rnite  Element  Method;  Phase  II 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 
Constnidion,  and  IMntsnanoe  of  Structures 

User's  Guide:  Conner  Program  for  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Investigation  and  Design  of  U-Franw  Structures  Using 
Program  CUFRBC 

Volume  A;  Program  Criteria  and  Documentation 
Volume  B:  User's  Guide  for  Basins 
Vofome  C:  User's  Guide  for  Channels 

User's  Guide:  Corr^er  Program  for  Two-Dimensional  Analysis 
of  U-Frame  or  W-Frame  Structures  (CWFRAM) 

User's  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPQC)  Preprocessor  to  CPGA  Program 

Application  of  Finite  Element,  Grid  Gerteration,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 
Groundwater  Modeling 

User's  Guide;  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Roare's  Momerri  Reduction 

User's  Guide  for  Concrete  Strength  investigation  and  Design 
(CASTR)  in  Accordance  with  ACI 31 8-89 

Rinite  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

Introduction  to  the  Computation  of  Response  Spectrum  for 
Earthquake  Loading 

Concept  Design  Example,  Computer  Aided  Structural 

Modeling  (CASM) 

Report  1 :  ^eme  A 
Report  2:  Scheme  6 
Report  3:  Scheme  C 

User's  Guide:  Computer-Aided  Structural  Modeling 

(CASM)  -  Version  3.00 

Tutorial  Guide:  Corr^xJter-Aided  Stmctural  Modeling 

(CASM)  -  Version  3.00 


Date 

Aug  1989 

Aug  1989 
Sep  1989 
Feb  1990 
May  1990 

Sep  1990 
Jun  1990 
Sep  1990 

Oct  1991 

Mar  1992 
May  1992 
Jun  1992 

Jun  1992 
Jun  1992 
Jun  1992 

Apr  1992 
Apr  1992 
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(Conduded) 

Title 


Date 


Contract  Report  m.-92-1 
Technical  Report  ITL-92-7 
Contract  Report  rTL-92-2 

Contract  Report  ITL-92-3 

Instruction  Report  GL-87-1 

Technical  Report  ITL-92-1 1 
Technical  Report  ITL-92-1 2 


Instruction  Report  GL-87-1 

Technical  Report  ITL-93-1 
Technical  Report  ITL-93-2 

Technical  Report  ITL-93-3 

Instruction  Report  ITL-93-3 
Instruction  Report  ITL-93-4 
Technical  Report  ITL-94-2 

Instruction  Report  ITL-94-1 

Instruction  Report  rTL-94-2 

Technical  Report  iTL-94-4 


Optimization  of  Steel  Pile  Foundations  Using  Optimality  Criteria 

Refined  Stress  Analysis  of  Melvin  Price  Locks  and  Dam 

Knowledge-Based  Expert  System  for  Selection  and  Design 
of  Retaining  Structures 

Evaluation  of  Thermal  and  Incremental  Construction  Effects 
for  Monoliths  AL-3  and  AL-5  of  the  Melvin  Price  Locks 
and  Dam 

User's  Guide:  UTEXAS3  Slope-Stability  Package;  Volume  IV, 

User's  Manual 


The  Seismic  Design  of  Waterfront  Retaining  Structures 


Computer-Aided,  Retd- Verified  Structural  Evaluation 

Report  1 :  Development  of  Computer  Modeling  T echniques 
for  MKsr  Lock  Gates 

Report  2:  Reid  Test  and  Analysis  Conelation  at  John  Hollis 
Bankhead  Lock  and  Dam 

Report  3:  RekJ  Test  and  Analysis  Correlation  of  a  Vertically 
Framed  Miter  Gate  at  Emsworth  Lock  and  Dam 


User's  Guide:  UTEXAS3  Slope-Stability  Package;  Volume  III, 
Example  Problems 


Theoretical  Manual  for  Analysis  of  Arch  Dams 

Steel  Structures  for  Civil  Works,  General  Considerations 
for  Design  and  Rehefoilitation 

Soil-Structure  Interaction  Study  of  Red  River  Lock  and  Dam 
No.  1  Subjected  to  Sediment  Loading 

User's  Manual — ADAP,  Graphics-Based  Dam  Analysis  Program 

Load  and  Resistance  Factor  Design  for  Steel  Miter  Gates 


User's  Guide  for  the  Incremental  Construction,  Soil-Structure  Interaction 
Program  SOILSTRUCT  with  Far-Reld  Boundary  Elements 

Tutorial  Guide:  Computer-Aided  Structural  Modeling  (CASM); 

Version  5.00 


User's  Guide:  Computer-Aided  Structural  Modeling  (CASM); 
Version  5.00 


Dynamics  of  Intake  Towers  and  Other  MDOF  Structures  Under 
Earthquake  Loads:  A  Computer-Aided  Approach 


Jun  1992 
Sep  1992 
Sep  1992 

Sep  1992 

Nov  1992 

Nov  1992 

Nov  1992 
Dec  1992 
Dec  1993 

Dec  1992 

Jul  1993 
Aug  1993 

Sep  1993 

Aug  1993 
Oct  1993 
Mar  1994 

Apr  1994 

Apr  1994 

Jul  1994 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the  originator. 


